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FOREWORD 


This  report  contains  the  results  of  a dissertation  which  was 
prepared  for  the  Aeronautical  and  Astronomical  Engineering  Depart- 
ment of  Ohio  State  University.  The  effort  was  supported  by  the 
Laser  Aerodynamics  Branch,  the  Components  Branch,  and  the  Interim 
Compressor  Research  Facility  Group  of  the  Air  Force  Aero  Propulsion 
Laboratory,  Air  Force  Systems  Command,  Wright-Patterson  AFB,  Ohio, 
under  Project  3012,  Task  14,  and  Work  Unit  0£.  The  effects  of  inte- 
gral scale,  micro  scale,  and  intensity  cf  turbulence  on  a laser  beam 
which  has  propagated  through  the  turbulent  medium,  as  measured  in  the 
focal  plane  of  the  laser  beam,  are  reported.  The  Geometric  ranges  of 
Turbulence  and  laser  beam  parameters  were  investigated  which  are  typi- 
cal of  laser  cavities.  The  effort  was  conducted  by  Richard  8.  Rivir 
during  the  period  of  March  1973  to  August  1976. 
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CHAPTER  I 


INTRODUCTION 

Since  the  development  of  the  first  flow  cooled  CC^  laser  in  1960, 
there  has  been  a steady  progression  in  the  development  of  flow  cooled 
lasers.  Pulsed  and  continuous  lasers  have  evolved  which  are  well  under- 
stood from  a physical  standpoint.  The  useful  applications  of  the  laser 
as  a precise  diagnostic  tool  in  flow  field  measurements,  as  a manufac- 
turing tool  for  drilling,  cutting,  and  welding  various  materials,  and  as 
a communications  or  power  link,  has  necessitated  higher  power  densities 
and  improved  beam  quality.  Increases  in  the  average  level  of  power  has 
increased  the  emphasis  on  the  use  of  a flowing  medium  to  separate  the 
high  entropy  lasing  medium  from  the  lasing  cavity.  The  continual  re- 
placement and  the  transport  of  the  lasing  medium  through  the  lasing 
cavity  introduces  fluctuations  in  the  physical  properties  of  the  lasing 
medium.  These  fluctuations  in  time  and  spatial  dimensions  of  the  physi- 
cal properties  of  the  fluid  can  be  of  the  same  general  magnitude  as  the 
optical  parameters  (time  and  spatial  dimensions)  of  the  laser  beam. 

When  optical  and  fluid  characteristic  times  and  dimensions  are  similar 
in  magnitude,  a coupling  may  occur  between  the  random  variations  in  the 
fluid  properties  and  the  electromagnetic  propagation  of  a wave  through 
the  medium. 
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Analytical  descriptions  of  the  coupling  of  the  electromagnetic 
field  and  the  fluctuating  fluid  properties  results  in  equations  which 
have  the  randomly  variable  fluid  properties  as  a coefficient  of  the 
sought  wave  function  solution.  An  exact  analytical  solution  to  this 
problem  is  not  known.  Analytic  approximations  which  make  use  of  vary- 
ing degrees  of  smallness  of  the  fluid,  transport,  physical,  and  optical 
properties  allow  qualitative  predictions  of  the  effect  of  the  randomly 
fluctuating  medium  on  the  electromagnetic  wave  propagation  through  the 
medium.  Current  analytical  models  developed  to  predict  the  degradation 
in  beam  quality,  as  a result  of  propagation  through  a turbulent  flow 
field,  yield  results  which  only  agree  within  an  order  of  magnitude  with 
the  observable  results. 

The  qualitative  effect  of  random  fluctuations  or  turbulence  on 
beam  spreading  is  well  demonstrated.  Turbulence  causes  scattering  of 
photons  as  a result  of  spatial  and  temporal  changes  in  the  index  of 
refraction  of  the  medium  with  a resultant  decrease  in  energy  per  unit 
area. 

The  increase  in  total  average  power  of  lasers  is  limited  basically 
by  economic  considerations  and  physical  space  limitations.  However, 
the  improvement  in  laser  beam  quality  and  average  power  density  (up  to 
atmospheric  breakdown)  is  dependent  in  part  on  a better  understanding 
of  the  interaction  between  the  fluctuations  of  the  physical  transport 
properties  of  the  fluid  through  which  the  beam  is  propagating  and  the 
physical  quantities  which  characterize  the  laser  beam.  The  degradation 
in  beam  quality  in  flow  cooled  lasers  is  within  the  accuracy  of  the 
theoretical  predictions  of  focal  spot  size  in  the  presence  of  turbulence. 


A better  correlation  of  theoretical  and  experimental  knowledge  of  the 
effects  of  turbulence  on  beam  quality  is  required  to  assess  the  possi- 
bility of  further  improvements  in  beam  quality  for  flow  cooled  lasers. 

The  fluctuations  of  physical  properties  of  the  lasing  fluid  or  a 
fluid  through  which  a beam  may  propagate  are  caused  by  the  interaction 
of  the  flow  with  an  assortment  of  physical  boundaries,  the  mixing  of 
dissimilar  component  gases  of  the  lasing  medium,  and  chemical  reactions. 

Examples  of  disturbance  generators  include:  nozzle  arrays,  mirrors, 

input  windows,  extraction  windows,  electrodes,  chemical  reaction  sources, 
thermal  sources,  wall  boundary  layers,  and  even  vortex  generators  which 
are  used  to  stabilize  some  electrical  discharge  lasers.  Specific  exam- 
ples will  be  illustrated  in  Chapter  II. 

The  number  of  theoretical  solutions  available  in  the  literature  for 
the  propagation  of  an  electromagnetic  wave  through  a randomly  varying 
medium  are  quite  numerous.  The  reason  for  many  theoretical  solutions  is 
due  to  the  lack  of  an  exact  solution  to  the  wave  equation.  The  wave 
equation  will  next  be  examined  and  the  solutions  grouped  into  generic 
types  as  they  are  found  in  the  literature.  This  section  will  then  be 
concluded  with  a historical  listing  of  some  of  the  publications  which 
deal  specifically  with  predicting  beam  spread  and  the  measurement  of 
beam  spread  in  the  presence  of  turbulence. 

The  intensity  amplitude,  intensity  distribution,  rms  level  of  the 
fluctuations  of  intensity,  and  the  spectrum  of  the  fluctuations  in  the 
presence  of  turbulence  are  quantities  of  interest  for  several  applica- 
tions and  for  which  appropriate  theoretical  predictions  are  possible. 

Each  of  these  quantities  is  a direct  measure  of  the  effects  of 

a 

3 

_J1 


turbulence.  However,  a clear  and  exact  relationship  between  the  causes 
and  the  observable  effects  does  not  exist. 

The  lack  of  an  exact  relationship  is  due  to  the  way  in  which  the 
random  fluctuations  of  the  medium  enter  the  propagation  equations.  The 
general  form  of  the  differential  equation  for  an  electromagnetic  wave 
propagating  in  a random  medium  is  obtained  by  eliminating  the  magnetic 
field  from  Maxwell's  equations  which  eventually  results  in: 

v E + k2  n2(r,t)  E + 2v  (E  • vln  n(r,t))=  0 (1) 

The  last  term  in  equation  (1)  represents  the  effects  of  depolariza- 
tion. Strohbehm  and  Clifford  have  calculated  that  this  contribution  to 
the  field  of  the  propagating  wave  is  160  db  less  than  the  incident  wave. 
Saleh  attempted  to  measure  this  component  with  a sensitivity  of  -45  db 
and  was  unsuccessful  (Lawrence  and  Strohbehm,  1970).  There  have  been  no 
attempts  to  estimate  the  contribution  to  the  depolarization  term  by  scat- 
tering. All  previous  analyses  reviewed  seek  a solution  for  the  electric 
field  which  neglected  the  depolarization  term.  Dropping  the  last  term 
results  in  the  scalar  wave  equation,  which  represents  the  simplest  theo- 
retical description  possible: 

v E + k2  n2(r,t)  E = 0 (2) 

The  random  fluctuations  are  seen  to  enter  the  above  equation  as  a 
coefficient  of  the  sought  solution  for  the  field  E.  This  is  the  source 
of  the  mathematical  difficulties,  since  there  is  no  known  exact  solution 
for  this  form  of  equation.  Various  approximations  and  degrees  of  small- 
ness are  required  along  with  their  implicit  ranges  of  validity  in  order 
to  obtain  even  approximate  solutions.  The  solutions  which  are  found  in 
the  literature  are  generally  characterized  by  one  of  the  following 
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methods,  or  range  of  geometric  parameters:  (1)  geometric  optics; 

(2)  method  of  small  perturbations;  (3)  method  of  smooth  perturbations 
(Rytov);  (4)  modulation  transfer  function  or  mutual  coherence  function; 
and  (5)  multiple  scattering.  Experimentally  usable  results  are  obtained 
from  the  first  four  methods.  All  of  the  above  methods  are  adequately 
documented  in  the  literature. 

The  first  three  methods  are  detailed  in  the  Russian  works  by 
Chernov  and  Tatarski,  which  were  translated  to  English  in  1960  and  1961 
respectively.  Many  of  the  papers  which  appear  in  the  literature  for  the 
following  eight  years  deal  with  modifications,  variations,  and  attempts 
to  confirm  observations  resulting  from  these  publications.  Huffnagal 
and  Stanley  introduced  the  fourth  method,  the  modulation  transfer 
function  and  an  average  coherence  function  in  an  analytical  solution 
published  in  1964.  In  1965,  Collins  and  Wortendyke  published  experi- 
mental scattering  measurements  for  a finite  plane  wave,  propagating 
through  a wind  tunnel,  which  were  based  upon  one  of  Tatarski's  solu- 
tions. In  1966,  Schmeltzer  published  an  iterative  solution  based  on 
Tatarski's  Rytov  approximation  (Method  3)  for  a focused,  initially 
gaussian  beam,  propagating  through  a turbulent  atmosphere.  Schmeltzer’s 
results  were  given  in  integral  form.  In  1968,  Gebhardt  evaluated  the 
integral  solutions  given  by  Schmeltzer  for  focused  and  colliminated 
gaussian  beams.  After  1969,  nearly  all  analytical  solutions  involve 
Method  4.  Sutton,  in  1969,  published  a calculation  of  a numerical  solu- 
tion using  the  modulation  transfer  function,  which  included  the  case  in 
which  the  turbulence  scale  and  the  beam  diameter  are  of  the  same  order. 
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In  1970,  Lutomirski  also  published  a solution  for  a finite  beam  for  some 
restricted  path  lengths  based  upon  the  mutual  coherence  function. 

There  are  several  excellent  surveys  currently  available  in  the  lit- 
erature. Lawrence  and  Strohbehm  published  "A  Survey  of  Clear-Air  Propa- 
gation Effects  Relevant  to  Optical  Communications"  in  the  Proceedings  of 
the  IEEE  in  1970.  An  extensive  survey  of  Russian  efforts  may  be  found 
in  the  1971  Soviet  Physics  article,  "Status  of  the  Theory  of  Propagation 
of  Waves  in  a Randomly  Inhomogeneous  Medium,"  by  Barabanenkov,  Kravtsov, 
Rytov,  and  Tamarskii . 

Lutomirski  and  Yura  published  numerous  reports  separately  and  to- 
gether. Most  of  their  reports  are  summarized  and  referenced  in  Rand 
publication,  "Degradation  of  Laser  Systems  by  Atmospheric  Turbulence," 
which  was  published  in  1973.  The  majority  of  the  solutions  presented 
by  the  Rand  Report  also  evolve  from  the  mutual  coherence  function.  Al- 
though this  effort  will  not  deal  with  scintillation  effects,  a very  good 
scintillation  review  is  given  by  Meyer-Arendt  and  Emmanuel  in  the  1965 
NBS  Technical  Report  225,  "Optical  Scintillation:  A Survey  of  the 

Literature. " 
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CHAPTER  II 

DEFINITION  OF  THE  PROBLEM 

In  this  chapter  some  examples  of  laser  cavity  flows  will  be  illus- 
trated to  show  the  types  of  turbulence  conditions  which  may  be  encoun- 
tered in  laser  cavity  flows.  The  lack  of  adequate  experimental  measure- 
ments to  confirm  existing  analytic  solutions  will  be  illustrated  as  the 
problem  to  be  investigated  in  laser  cavity-like  flows.  Next,  the  ex- 
perimental conditions  will  be  examined  and  selected  so  that  beam  spreads 
of  the  order  of  three  times  that  of  a diffraction  limited  beam  may  be 
produced  for  geometric  conditions  such  as  are  encountered  in  laser 
cavities.  Finally,  the  scope  or  range  of  experimental  conditions  to  be 
investigated  will  be  stated. 

Examples  of  Laser  Cavity  Flows 

Examples  of  simulated  laser  cavity  flows  are  shown  in  Figures  1 
and  2.  The  first  figure  represents  the  proposed  injection  of  the  lasing 
medium  for  a pulsed  electrical  discharge  laser  (EDL).  The  first  figure, 
a schlieren  photograph,  shows  the  injection  of  the  lasing  medium  from  a 
series  of  parallel  manifolds  which  contain  a large  number  of  small  jets. 
The  pressure  ratio  across  the  jet  orifices  is  greater  than  critical. 

The  shear  layers  surrounding  the  jets  are  seen  to  remain  well  defined 
for  approximately  100  jet  diameters  downstream.  This  ordered  level  of 
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Figure  1.  Simulated  Cold  Flow  of  a Proposed 
Injection  Technique  for  an  EDL 
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turbulence  then  decays  into  a region  of  isotropic  turbulence  in  the 
upper  portion  of  the  figure  where  the  scale  size  of  the  turbulence  is 
observed  to  be  approximately  the  same  in  all  directions.  A further 
problem  encountered  with  large  numbers  of  jets,  which  is  not  specifi- 
cally illustrated  by  Figure  1,  is  the  interaction  of  individual  jets 
to  produce  both  unsteady  motion  and  various  degrees  of  coalescence  by 
the  interaction  of  individual  jets,  thereby  increasing  the  integral 
scale  of  turbulence  which  is  observed.  The  lasing  path  for  this  par- 
ticular device  would  be  in  the  plane  of  the  photograph  and  normal  to 
flow.  The  pressure  of  the  lasing  medium  in  an  electrical  discharge 
laser  would  typically  be  on  the  order  of  an  atmosphere  and  the  mean 
flow  would  be  at  a Mach  number  on  the  order  of  .1.  The  turbulence 
Reynolds  number  of  this  flow  would  be  on  the  order  of  7000.  For  injec- 
tion techniques  of  this  type,  the  lasing  path  would  include  regions  of 
both  isotropic  and  nonisotropic  turbulence.  The  turbulence  intensities 
in  the  nonisotropic  regions  of  supercritical  jets,  such  as  shown  in 
Figure  1,  are  on  the  order  of  50%;  as  the  region  approaches  isotropic 
turbulence,  the  turbulence  intensity  decreases  to  a few  percent. 

The  second  figure  shows  a multiple  element  supersonic  nozzle  array 
similar  to  that  encountered  in  gas  dynamic  (GDL)  and  chemical  lasers 
where  the  throat  diameters  are  on  the  order  of  a few  hundredths  to  a 
few  thousandths  of  a centimeter.  This  schlieren  photograph  shows  Mach 
waves  and  a shear  layer  leaving  each  nozzle.  The  combination  of  a low 
lasing  medium  pressure  and  the  very  small  characteristic  lengths  of  the 
throat  and  boundary  layer  results  in  a value  of  turbulence  Reynolds 
number  on  the  order  of  a few  hundred.  A comparison  of  Figures  1 and  2 
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shows  that  the  wakes  are  propagated  for  much  longer  distances  as  viscous 
and  inertial  forces  approach  one  another.  The  wakes  are  terminated  only 
by  the  formation  of  regions  of  even  greater  turbulence  intensity.  Simi- 
lar weaker  wakes  will  also  be  observed  at  the  intersection  of  shocks  of 
unequal  strength  or  turning  angle.  The  velocity  turbulence  intensity  in 
these  wakes  can  again  approach  50%. 

In  AIAA  Paper  73-626  (Director,  1973),  density  variations  (Ap/p) 
observed  in  simulated  cavity  flows  from  wakes  and  wall  contours  with 
interferograms  range  from  a few  tenths  of  a percent  to  approximately 
five  percent  for  flows  similar  to  that  shown  in  Figure  2.  The  displace- 
ment boundary  layer  thickness,  or  one  half  wake  thickness,  was  calcu- 
lated to  be  on  the  order  of  .05  centimeter  in  cold  flow.  In  AIAA  Paper 
73-641  (Hyde  and  Hosack,  1973),  the  thermal  boundary  layer  was  calcu- 
lated for  a chemical  laser  flow  nozzle.  The  energy  boundary  layer  was 
observed  to  be  four  times  the  thickness  of  the  velocity  boundary  layer, 
on  the  order  of  .76  to  1 centimeter,  a very  significant  portion  of  the 
exit  area  of  the  nozzle,  even  with  nozzle  area  ratios  which  are  quite 
large.  The  thermal  wake  is  the  region  of  interest  in  terms  of  its 
influence  on  the  ability  of  turbulence  to  introduce  disturbances  and 
defocus  a beam  propagating  through  the  medium. 

In  the  electro  aerodynamic  laser  (EAL),  a series  of  half  cylinders 
are  arranged  in  a picket  fence  fashion  as  the  anode  support  structure 
at  the  entrance  of  the  lasing  cavity.  Each  of  these  cylinders  sheds 
vortices  or  streaks  periodically.  The  turbulence  Reynolds  number  of  the 
flow  in  the  electrical  discharge  laser  cavity  is  on  the  order  of  a few 
thousand  due  to  the  low  density  and  also  low  velocity.  Since  the 
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Reynolds  number  is  low,  the  shedding  structure  embedded  in  the  flow 
remains  for  long  distances  just  as  the  wakes  do  in  the  GDL  and  EAL. 

The  low  pressure  lasing  mediums  in  cavities  such  as  the  GDL,  EAL, 
and  chemical  lasers  may  also  require  aerodynamic  windows.  These  windows 
may  involve  shocks  and  the  intersection  of  complex  patterns  of  shock 
waves.  At  the  intersection  of  complex  or  unequal  strength  shock  waves, 
complex  regions  of  turbulence  also  exist. 

The  turbulence  Reynolds  numbers  are  observed  to  be  very  small  in 
laser  cavities  compared  to  the  turbulence  Reynolds  numbers  encountered 
in  atmospheric  propagation.  The  Reynolds  number  in  atmospheric  propaga- 
tion problems  are  large  due  to  the  scale  of  the  disturbances.  In  cavi- 
ties the  beam  diameter  to  turbulence  scale  can  range  from  much  less  than 
1 to  ^1 . Turbulence  intensities  may  be  quite  high  (^50%)  and  the  turbu- 
lence conditions  range  from  nonisotropic  to  isotropic.  The  path  length 
through  the  lasing  medium  is  short  in  comparison  to  atmospheric  propaga- 
tion problems,  ranging  only  to  a few  tens  of  meters.  The  lasing  beam  is 
near  gaussian  or  enters  the  cavity  as  a gaussian  beam.  The  lasing  beam 
would  normally  be  a constant  diameter,  although  on  long  path  length 
devices,  some  refocusing  is  desirable.  If  the  laser  is  to  be  utilized 
at  its  focal  plane,  the  spread  of  the  beam  in  the  focal  plane,  due  to 
turbulence  in  the  lasing  medium,  is  of  interest. 

Examples  of  Correlation  of  Experiment  with  Theory  and  Comparison  of 
Similar  Theories 

In  the  literature,  the  majority  of  the  publications  on  the  propaga- 
tion of  an  electromagnetic  wave  through  a randomly  varying  medium  are  of 

I 


11 


a theoretical  nature.  This  experiment  will  be  concerned  with  only  the 
prediction  and  measurement  of  beam  spread.  The  1965  paper  by  Collins 
and  Wortendyke  reported  one  of  the  efforts  to  correlate  the  propagation 
of  finite  rectangular  plane  wave  with  one  of  Tatarski’s  solutions.  This 
effort  is  one  of  the  better  attempts  to  make  measurements  in  a control- 
led environment.  The  turbulence  properties  of  the  flow  field  were  not 
measured,  but  calculated  from  the  results  of  a similar  experiment  in 
another  facility.  The  claim  of  Collins  and  Wortendyke's  work  was  order 
of  magnitude  agreement  with  theory,  with  factors  typically  ranging  from 
two  to  seven  times  the  theoretical  value  of  the  power  scattered.  Exper- 
imental results  are  presented  of  power  scattered  or  spread  as  a function 
of  mean  temperature  rise,  angle  of  scattering,  and  scattering  volume. 

In  general,  the  slopes  of  the  curves  in  addition  to  differing  in  abso- 
lute value  from  the  analytical  predictions,  they  also  differ  in  the 
relative  amount  of  parametric  change  in  the  experimental  curves  from  the 
analytical  predictions.  Further  examination  of  the  experimental  results 
also  shows  that  two  of  the  three  experimental  curves  for  the  volume 
dependence  of  beam  spread  have  the  opposite  slope  of  the  theoretical 
curves.  There  exists  in  the  literature  very  weak  experimental  documen- 
tation of  even  the  trends  of  the  simplest  analytical  models.  There  are 
virtually  no  thorough  experimental  investigations  of  the  limits  of  appli- 
cability of  the  theoretical  models.  As  a result,  when  examining  much  of 
the  experimental  work,  it  is  not  clear  whether  the  theoretical  models 
are  inadequate  or  whether  the  lack  of  accurate  turbulence  measurements 
precludes  an  accurate  prediction. 
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A report  which  should  be  included  in  the  literature  surveys  is 
the  1971  Riverside  Research  Institute's,  "How  Well  Can  Laser  Radiation 
be  Focused  Through  the  Atmosphere,"  by  Greenbaum.  Greenbaum  points 
out  the  significantly  different  theoretical  predictions  of  spot  size 
obtained  by  Gebhardt  and  Collins  (Gebhardt  and  Collins,  1969,  a paper 
based  on  Gebhardt's  1968  paper)  and  that  of  Freid  (Freid,  1966). 
Greenbaum's  figure  comparing  the  two  theories  is  shown  in  Figure  3. 
Greenbaum' s Reference  3 curves  refer  to  Gebhardt  and  Collins'  work  and 
his  Reference  4 refers  to  Freid's  work.  The  possibility  of  more  than 
order  of  magnitude  differences  between  similar  theories  is  clearly 
indicated  by  Greenbaum's  comparison. 


Figure  3.  Greenbaum's  Comparison  of  Required  Initial  Beam 
Radii  Given  Maximum  Allowable  Focal  Spot  Diameter 
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It  is  apparent  that  the  experimental  data  currently  available  are 
inadequate  so  that  neither  the  trends  nor  the  quantitative  predictions 
of  the  numerous  analytical  approaches  which  are  available  can  be  judged. 
The  experimental  information  which  is  missing  most  often  is  an  accurate 
description  of  the  turbulent  flow  field.  This  is  in  part  due  to  tne 
fact  that  design  and  control  of  experimental  parameters  to  cover  all 
ranges  and  explore  the  boundaries  of  even  one  theory's  variables  is  a 
sizable  undertaking. 

Since  the  initial  interest  in  this  problem  arose  from  the  desire 
to  accurately  predict  what  portion  of  the  beam  spread  of  fast  flow 
lasers  might  be  due  to  turbulence,  within  the  lasing  cavity,  the  selec- 
tion of  the  experimental  parameters  will  be  based  upon  those  typical  of 
laser  cavities  rather  than  the  values  which  are  necessary  to  prove  or 
disprove  any  one  particular  theory. 

Statement  of  the  Problem 

The  objective  of  this  effort  is  to  obtain  adequate  experimental 
data,  both  in  definition  of  the  turbulent  flow  field  and  in  beam  spread, 
to  determine  if  the  discrepancies  between  experiment  and  theory  for  the 
prediction  of  the  focal  spot  size  are  due  to  inadequate  turbulence 
measurements  or  inadequate  theories.  These  measurements  will  be  accom- 
plished in  a controlled  turbulent  flow  field  under  geometrical  con- 
straints which  are  reasonably  typical  of  fast  flow  laser  cavities. 

Selection  of  Experimental  Parameters  and  Their  Rang e s 

The  parameters  which  appear  in  the  theoretical  expressions  (Chernov, 
1960;  Tatarski,  1961;  Freid,  1966;  Gebhardt,  1968;  Gebhardt  and  Collins, 
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1969;  Sutton,  1969;  Lutomirski,  1970;  Lutomirski  et  al  , 1973)  for  beam 
spread  include  the  fluctuation  in  the  index  of  refraction  <An>,  the 
laser  beam  diameter  to  wave  length  ratio,  the  laser  beam  diameter  to 
integral  scale  of  turbulence  ratio,  the  path  length,  and  under  some  con- 
ditions, the  ratio  of  laser  beam  diameter  to  the  microscale  of  turbu- 
lence (Lutomirski  et  al , 1973).  For  integral  formulation,  the  ratios  to 
the  beam  diameter  to  wave  length  or  turbulence  scale  are  not  always 
readily  apparent  (Lutomirski,  197u). 

Initially,  consideration  was  given  to  making  turbulence  measure- 
ments in  a fast  flow  CO2  laser  cavity.  The  expected  range  of  experi- 
mental conditions  is  then  rather  limited  to  a specific  laser's  limited 
range  of  operation.  Beam  quality  can  also  be  coupled  with  the  particu- 
lar operating  point. 

While  cold  flow  turbulence  measurements  can  be  obtained  in  actual 
laser  cavities  with  conventional  techniques,  hot  flow  turbulence  measure- 
ments or  measurements  with  excitation  and  lasing  of  the  flowing  medium 
would  be  exceedingly  difficult  to  obtain.  Since  analytical  predictions 
and  experimental  measurements  currently  differ  by  an  order  of  magnitude, 
the  confirmation  of  the  predicted  analytical  trends  are  the  first  issues 
to  be  resolved.  Selection  of  a very  simple  and  controllable  turbulent 
flow  field  is  in  order  to  eliminate  the  uncertainties  associated  with 
unconventional  measurement  techniques  which  would  be  required  in  a laser 
cavity. 

The  investigation  of  focal  spot  size  for  laser  wave  lengths,  such 
as  CO2,  requires  additional  analytic  considerations  since  the  wave 
length  and  the  scale  of  turbulence  are  of  the  same  order  of  magnitude. 


15 


Once  again,  it  must  be  considered  that  even  when  the  wave  length  is  much 
smaller  than  the  turbulence  scale  size  (in  the  visible  spectrum  of  wave 
lengths),  the  analytic  prediction  and  experimental  measurements  are  at 
best  only  within  an  order  of  magnitude  of  each  other,  and  even  then 
trends  are  not  necessarily  consistent  with  theoretical  predictions  as 
was  shown  previously.  The  use  of  a CO2  wave  length  and  the  associated 
complications  with  windows,  stability,  alignment,  and  safety  requirements 
are  not  necessary  to  resolve  the  fundamental  issues.  Two  wave  lengths. 

He  Ne,  532. 8nm,  and  He  Cd,  441. 6nm,  were  selected  for  this  investigation. 

It  was  previously  indicated  in  Examples  of  Laser  Cavity  Flows  that 
the  ratio  of  laser  beam  diameters  to  turbulence  scale  size  ranges  from 
:<1  to  'vl.  Figure  4 shows  the  range  of  beam  diameters  which  are  possi- 
ble with  a conventional  Spectra  Physics  col  1 iminator  and  beam  expanding 
telescope  for  He  Ne  and  He  Cd.  Also  shown  in  Figure  4 are  the  compara- 
ble curves  for  HF,  DF,  CO,  and  CO2  lasing  mediums  for  the  range  of  beam 
diameter  to  wave  length  which  can  be  obtained  from  the  He  Ne  and  He  Cd 
beams  with  the  conventional  Spectra  Physics  beam  expanding  telescope. 
Apertures  and  lenses  were  selected  for  the  Spectra  Physics  beam  expander 
:o  give  nominally  50mm,  26mm,  and  8mm  He  Ne  beams. 

The  scale  of  turbulence  is  controlled  by  the  physical  dimensions  of 
the  turbulence  generator,  the  diameter  of  a rod,  the  diameter  of  the 
perforation,  or  the  mesh  size.  The  turbulence  Reynolds  numbers  typical 
of  laser  cavities  were  indicated  (in  Examples  of  Laser  Cavity  Flows)  to 
be  quite  low.  The  integral  scale  of  Reynolds  number  will  later  be  shown 
to  be  related  to  the  microscale  by  the  relation: 


(3) 
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where  Ag  is  a constant  which  ranges  between  4 and  30.  Relation  (3)  is 
displayed  in  Figure  5.  The  range  of  microscales  to  integral  scales, 
which  would  be  expected  to  be  found  in  the  different  types  of  laser  cav- 
ities, has  been  indicated  in  Figure  5.  The  Reynolds  numbers  are  seen  to 
run  from  the  critical  range  to  20,000.  The  EAL  has  integral  Reynolds 
numbers  in  the  range  from  300  to  8,000  due  to  the  low  velocities  and  cav- 
ity densities.  Typical  pulsed  laser  cavities  exhibit  integral  Reynolds 
numbers  on  the  order  of  10,000.  Even  the  wakes  from  GDL  nozzles  are  seen 
to  be  in  the  range  of  200  to  20,000.  The  size  of  the  wake  shed  by  a 
nozzle  is  proportional  to  the  boundary  layer  thickness.  Scaling  of  the 
boundary  layer  thickness  is  Re  to  Re  , depending  on  whether  the 
flow  is  laminar  or  turbulent.  Since  the  boundary  layer  in  the  nozzles 
is  quite  small,  the  Reynolds  number  is  likewise,  quite  small.  Simula- 
tion of  turbulence  Reynolds  numbers  in  the  range  of  100  to  10^  is  of 
interest  in  terms  of  scaling  the  effects  of  the  integral  and  the  micro- 
scales on  beam  quality  for  the  GDL.  Pulsed  lasers  generally  operate  at 
higher  pressures  and  turbulence  Reynolds  numbers  are  on  the  order  of 
9,000  to  17,000.  The  range  of  turbulence  Reynolds  numbers,  which  are  of 
interest  fo^  cavity  propagation  problems,  ranges  from  100  to  20,000. 

The  range  of  turbulence  Reynolds  numbers  which  are  feasible  in  the 
AFIT  23  centimeter  wind  tunnel  are  indicated  in  Figure  6.  Figure  6 shows 
lines  of  constant  velocity,  on  which  the  tunnel  is  capable  of  operating. 
The  tunnel  will  control  and  repeat  tunnel  velocity  settings  from  .61 
meters  per  second  to  16.8  meters  per  second  with  an  empty  test  section. 
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Figure  5.  Laser  Cavity  Scales  of  Turbulence 


Turbulence  generators  with  characteri stic  diameters  of  .045,  .64,  and 
1.27  centimeters  were  chosen  as  the  dimensions  of  the  turbulence  genera- 
tors for  this  work.  The  scales  associated  with  these  generators  are 

indicated  by  the  horizontal  lines  across  the  range  of  available  tunnel 
velocities  in  Figure  6. 

The  criteria  which  are  used  to  justify  the  neglect  of  diffraction 
effects  state  that  the  microscale  iQ  must  be  greater  than  the  square 
root  of  the  product  of  the  wave  length  and  the  path  length,  Alq.  Fig- 
ure 7 displays  the  range  of  AL0  in  typical  laser  cavities  for  the 
lasing  mediums  of  He  Cd,  He  Ne,  HF,  CO,  and  CO2.  The  /AL0  ranges  from 
.5  centimeter  to  13  centimeters  for  HF  and  from  1 to  30  centimeters  for 

CO2.  The  figure  also  indicates  that  a microscale  greater  than  .5  centi- 

meter is  required  in  the  AFIT  wind  tunnel  in  order  to  neglect  the  dif- 
fraction effects  for  helium  neon  and  helium  cadmium  lasers.  Although 
the  results  from  this  effort  are  not  intended  to  include  any  considera- 
tion of  diffractive  effects,  the  results  may  nevertheless  have  limited 
direct  implications  for  fast  flow  laser  applications. 

The  path  length  and  the  fluctuation  in  index  of  refraction  remain 
as  experimental  parameters  to  be  specified.  A single  pass  path  length 
of  23  centimeters  was  available  in  the  AFIT  tunnel;  however,  multiple 
passes  could  be  achieved.  Only  a single  pass  will  be  reported  in  this 
work,  although  limited  data  were  also  collected  for  three  and  five 
passes.  The  fluctuations  in  index  of  refraction  are  introduced  by 
heated  wires  such  that  beam  spreads  up  to  three  times  a diffraction 
limited  beam  could  be  achieved. 


w 


20 


|S«BKK 


Ktti 


U 

CD 


+->  O') 
cr>  d) 


itions  for  Geometric  Optics 
Diffractive  Effects 


Scope  of  the  Investigation 


By  varying  optical  and  fluidic  lengths,  temperature,  and  temporal 
scales,  the  range  of  flow  parameters  applicable  to  fast  flow  laser  cavi- 
ties will  be  considered  in  order  to  determine  whether  the  analytical 
predictions  are  functionally  correct,  accurate,  and  to  improve  the  base 
of  experimental  information  from  which  the  predictions  of  beam  spot  size 
may  be  formulated  for  laser  cavities.  The  experimental  investigation 
will  emphasize  low  values  of  turbulence  Reynolds  numbers  where  there  is 
very  limited  experimental  data.  The  fluctuations  in  index  of  refraction 
will  range  through  strong  turbulence  conditions. 

The  range  of  parameters  which  will  be  investigated  are  as  follows: 

1.  Mean  fluid  velocity  ranged  from  zero  to  20  meters  per  second 

2.  Turbulence  Reynolds  numbers  to  16,000 

3.  Velocity  turbulence  intensities  to  40% 

4.  Temperature  turbulence  intensities  to  .6% 

5.  Index  of  refraction  fluctuations  to  3.7  x 10"^ 

6.  Propagation  path  length  - 23  cm 

7.  Nominal  laser  beam  diameters  - 50mm,  30nm,  14mm,  and  1.5mm 

8.  Wave  lengths  of  the  laser  beams  of  632. 8nm  and  441. 6nm 

9.  Beam  diameter/turbulence  scale  - D/A  <<1  to  D/A  ^1 

10.  Attenuation  coefficient  (to  be  defined  in  Chapter  IV)  - 
OiLq  <<1 , «Lq  ^<1 , aL0  >1 
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CHAPTER  III 

APPARATUS  AND  INSTRUMENTATION 

Schematic  of  the  Experimental  Apparatus 

The  experimental  set-up  will  be  traced  from  the  open  circuit  wind 
tunnel  through  the  instrumentation  employed  for  data  acquisition,  in- 
cluding the  on-line  real  time  computer  used  for  calculation  of  the  cor- 
relation function.  Photographs  identifying  the  major  components  of  the 
apparatus  and  instrumentation  are  shown  in  Figures  8 and  9.  The  laser 
is  mounted  either  normal  or  parallel  to  the  longitude  tunnel  centerline 
and  then  turned  through  90°  on  the  tunnel  horizontal  centerline.  Figure 
8 shows  the  50  milliwatt  He  Cd  laser  in  position.  The  laser  beam  is 
expanded  and  spatially  filtered.  The  output  beam  is  then  passed  through 
a final  aperture  of  the  desired  diameter.  The  laser  beam  may  be  posi- 
tioned +10  centimeters  along  the  tunnel  centerline  by  moving  either  the 
laser  or  the  90°  turning  mirror  on  an  optical  bench  parallel  to  the 
tunnel  axis.  A schematic  of  the  experiment  and  the  optical  path  is 
shown  in  Figure  10.  A template  was  used  to  relocate  and  position  the 
beam  to  the  same  geometric  location  in  the  tunnel.  The  expanded  laser 
beam  enters  and  leaves  the  turbulent  medium  through  interferometer 
quality  windows  which  are  parallel  and  flush  with  the  tunnel  walls.  The 
turbulence  is  controlled  and  generated  at  a fixed  location  upstream  of 
the  tunnel  windows.  The  beam,  upon  emerging  from  the  tunnel,  enters  a j 
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Figure  10.  Schematic 


light  pipe  in  which  it  makes  six  traverses  in  order  to  reach  the  far 
field  for  the  largest  beam  diameter  used.  Although  near  field  measure- 
ments or  data  are  also  of  interest  in  laser  cavities,  only  far  field 
measurements  will  be  reported  in  this  work.  The  rationale  being  to  exam- 
ine the  effect  at  a target  due  to  cavity  disturbances.  The  beam  for  all 
measurements  to  be  reported  in  this  work  was  focused  at  45.7  meters.  At 
the  focal  plane,  an  aperture  is  traversed  across  the  plane  and  the  trans- 
mitted power  measured  with  a thermopile.  The  spatial  distribution  of 
the  energy  spectrum  and  autocorrelation  of  the  laser  beam  was  measured 
for  the  laser  beam  at  the  focal  plane.  The  turbulence  parameters  were 


obtained  by  double  parallel  hot  wires.  The  wires  were  operated  at 


different  overheat  ratios  to  obtain  mean  value,  turbulence  intensities, 
tunnel  profile,  autocorrelation,  and  energy  spectrums  of  velocity  and 
temperature. 

Wind  Tunnel 

An  open  circuit  wind  tunnel,  with  a 22.86  centimeter  by  22.86 
centimeter  test  cabin,  one  meter  in  length,  provided  the  controlled  flow 
field  for  the  research  (Figure  9).  The  room  which  cc"ta’r'.ed  the  wind 
tunnel  acted  as  the  closing  element  for  the  tunnel  and  provided  a heat 
sink  for  energy  dissipation  during  the  operation  of  heated  vortex  gener- 
ators. The  tunnel  has  a contraction  ratio  of  100  to  1.  With  a clean 
test  cabin,  the  tunnel  is  capable  of  velocities  of  18.34  meters  per 
second,  with  a turbulence  intensity  level  (u'/U)  of  10“3  and  a wall 
boundary  layer  0.025  meters  thick  at  the  mid-point  of  the  test  cabin. 

The  test  cabin  was  fitted  with  two  interferometer  quality  windows  0.203 
meters  in  diameter,  isolated  and  independently  suspended  from  the  tunnel 
by  a mounting  on  an  optical  bench.  The  windows  were  aligned  parallel  to 
the  flow  and  held  to  a divergence  angle  less  than  one  part  in  7000. 

Vortex  Generators 

Two  vortex  generator  systems  were  utilized  to  introduce  random 
fluctuations  into  the  flow  field.  First,  a perforated  plate  mounted 
normal  to  the  flow  upstream  of  the  mid-point  of  the  test  cabin  provided 
controlled  velocity  fluctuation.  Second,  an  array  of  parallel  nichrome 
wires  normal  to  the  flow  and  located  between  the  mid-point  of  the  test 
cabin  and  the  perforated  plate  provided  temperature  fluctuation. 
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Two  different  perforated  plates  were  used  to  provide  fluctuations 
in  the  velocity  field  and  variation  in  turbulence  scale.  Plate  A was 
0.317  centimeter  thick  with  a sixty  degree  array  of  .197  centimeter 
diameter  holes  with  a blockage  of  52.5  percent.  The  plate  was  located 
14  hole  diameters  upstream  from  the  test  cabin  mid-point.  The  other 
plate,  Plate  B,  was  0.317  centimeter  thick  with  a sixty  degree  array  of 
.098  centimeter  diameter  holes  and  a blockage  of  65.2  percent.  The 
plate  was  located  34.5  hole  diameters  upstream  from  the  test  cabin  mid- 
point. 

The  empirical  relationship  between  the  microscale  and  the  integral 
scale  was  given  previously  by  equation  (3)  as: 


VA 
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The  scale  size  of  turbulence  may  be  controlled  experimentally  by 
varying  the  mesh  size  or  the  characteristic  length  of  the  shedding  ele- 
ment and  the  free  stream  velocity  as  is  shown  by  the  equation  above  for 
the  microscale.  The  relationship  between  integral  and  microscale  may 
be  controlled  experimentally  as  shown  by  the  empirical  equation  (3)  by 
varying  the  Reynolds  number  of  the  flow.  The  dependence  of  the  integral 
scale  upon  the  free  stream  velocity  and  the  characteristic  length  of  a 
shedding  element  is  illustrated  by  the  Strouhal  relationship: 

f = Dp/U  (4) 

where  f is  the  frequency,  S the  strouhal  number,  and  Dp  the  diameter  of 
the  shedding  element. 

Experimental  control  of  temperature  fluctuations  is  required, 
since  the  primary  effect  of  random  fluctuations  on  the  laser  beam  being 


propagated  through  the  medium  results  from  flucutations  in  temperature 
which  cause  fl uctuation  in  the  index  of  refraction.  Temperature  fluc- 
tuations were  introduced  into  the  tunnel  by  a grid  of  electrically 
heated  ni chrome  wires.  The  wires  were  .046  centimeter  in  diameter,  on 
2.54  centimeter  centers,  beginning  2.54  centimeters  from  the  tunnel 
upper  wall.  The  power  input  into  the  wire  grid  ranged  from  .57  to  ^1.2 
kilowatts,  or  a mean  temperature  rise  of  the  flow  field  of  up  to  53°C. 
The  wire  grid  was  located  10.16  centimeters  downstream  from  the  perfo- 
rated plate  and  in  a plane  normal  to  the  windows.  Temperature  fluctua- 
tions, although  they  decay  slower  than  velocity  fluctuations,  show  sim- 
ilar relationships  to  those  of  velocity  fluctuations.  Similar  decay 
curves  for  both  velocity  and  temperature  have  been  measured  (Mills, 
Kistler,  Obrien  and  Corrsin,  1958,  p.  47,  Figure  17,  and  p.  50,  Figure 
19).  The  decay  rate  of  the  ratio  of  velocity  to  the  temperature  field 
approaches  .6  for  large  values  of  x/Dp.  The  ratio  of  integral  tempera- 
ture scale  to  the  integral  velocity  scale  is  greater  than  .9  for  temp- 
erature and  velocity  turbulence  generators  which  have  the  same  charac- 
teristic diameter  (Mills,  Kistler,  Obrien  and  Corrsin , 1958,  p.  66, 
Figure  30). 

Anemometer  System 

Mean  and  rms  velocity  and  temperature  profiles  of  the  tunnel  were 
obtained  with  a Thermo ‘Systems  (TSI)  Model  1050  anemometer,  Model  1057 
signal  conditioner,  and  Model  1052  polynominal  linearizer.  The  signal 
conditioner  was  used  primarily  as  an  amplifier  for  cases  of  very  small 
temperature  fluctuations.  The  linearizer  was  used  only  for  x probe 
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measurements  since  higher  accuracies  are  achievable  with  direct  read- 
ings, particularly  at  low  velocities  ('vl  m/sec).  The  measurement  of 
temperature  and  velocity  will  be  discussed  in  detail  in  Chapter  IV. 
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The  local  temperatures  and  velocities  were  obtained  with  parallel 
wires.  Thermo  Systems  Probe  Model  1244  T 1.5.  During  optical  measure- 
ment, the  wires  were  located  at  the  downstream  edge  of  the  beam  on  the 
horizontal  centerline  of  the  beam  and  at  the  lateral  center  of  the 
tunnel.  Tunnel  traverses  will  be  reported  with  the  wire  at  the  beam 

1 

centerline  locatiom— — _ 

Optical  Components 
Optical  Bench 

A three-ton  cast  iron  bed  plate,  two  by  three  meters,  isolated 
from  building  vibrations  by  2.54  centimeter  thick  pads  of  40  durometer 
rubber,  was  used  as  a table  to  support  optical  benches.  Two  optical 
benches  were  utilized.  Bench  C supported  the  tunnel  windows  and  the 
helium  neon  laser.  Bench  D supported  the  helium  cadmium  laser.  Bench 
D was  necessary  due  to  the  size  of  the  helium  cadmium  laser  and  the  lack 
of  room  normal  to  the  propagation  path.  Bench  E was  a two-ton  bed  plate 
(1  by  2.5  meters).  It  was  located  seven  meters  from  the  wind  tunnel  on 
a line  normal  to  the  tunnel  at  the  test  cabin  mid-section.  The  table 

: 

supported  the  common  mirror  used  to  fold  the  laser  beam  oack  and  forth 
across  the  room. 

i 

Laser  Sources 

A 15  milliwatt  Spectra  Physics  Model  124A  helium  neon  and  a 
50  milliwatt  Spectra  Physics  Model  185  helium  cadmium  laser  were  used 
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as  the  sources.  The  long  term  drift  for  the  sources  was  indicated  by 
the  manufacturer  to  be  less  than  5%  and  the  beam  amplitude  noise  and 
ripple  to  be  less  than  .3%  rms  from  120  Hertz  to  100  Kilohertz.  The 
output  beam  diameter  at  l/e:>  for  the  helium  neon  laser  is  1.1  milli- 
meters and  1.5  millimeters  for  the  helium  cadmium  laser.  The  output 
beams  of  these  lasers  were  expanded  with  the  Spectra  Physics  Model  336 
hi  332  collimating  lens,  expanding  lens,  and  spatial  filter  assembly. 
Nominal  beam  diameters  of  50  millimeters,  26  millimeters,  and  8 milli- 
meters in  diameter  were  obtained  for  the  helium  neon  and  helium  cadmium 
lasers.  Results  with  the  1.5  millimeter  helium  cadmium  beam  will  also 
? reported.  The  laser  beam  was  focused  by  the  Model  336  collimating 
lens  at  a distance  of  45.70  meters.  The  beam  remains  approximately 
onstant  in  cross  section  across  the  tunnel,  ^enter^exit  ^or  ^ 
illimeter  diameter  beam  is  1.0025.  The  50  millimeter  diameter  beam 
i filled  the  collimating  lens  such  that  the  lens  assembly  acts  as  a 
millimeter  aperture.  Additional  exit  apertures  were  used  for  the 
m,  i 1 er  diameter  beams  to  reduce  the  fringing  and  provide  a clean 
u.sian  input  beam  to  the  tunnel  test  section. 

Ihe  Lujht  Pijae 

Because  the  total  optical  path  was  200  times  the  turbulent 
I i urn  path,  a light  pipe  was  required  to  reduce  the  background  noise 
ceptable  levels.  The  laser  beam  path  length,  measured  from  the 
of  the  wind  tunnel  test  cabin  to  the  focal  plane,  ranged  from  6 to 
meters.  This  path  length  was  provided  by  folding  the  beam  inside 
light  pipe.  Three  light  pipes  were  utilized  during  the  course  of 
investigation.  Light  Pipe  A was  a 20  centimeter  diameter  aluminum 
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duct  seven  meters  long.  Due  to  the  limited  laser  beam  space  in  Pipe  A, 
a second  light  pipe,  Pipe  B,  was  built.  It  was  a .92  by  .92  meter 
rectangular  duct,  seven  meters  long,  fabricated  from  visqueen  wrapped 
on  a wood  framing.  In  order  to  enhance  the  signal  to  noise  ratio  of 
the  output  of  the  laser  beam  and  minimize  the  breathing  observed  with 
the  visqueen,  the  visqueen  light  pipe  was  replaced  by  a solid  structure 
fabricated  from  masonite  as  shown  in  Figure  7.  After  enclosure  of  the 
external  optical  path,  the  rms  value  of  laser  fluctuations  at  the  focal 
spot,  after  traversing  the  external  optical  path,  was  nominally  within 
the  laser  stability  specifications  of  .3%.  Both  the  visqueen  and 
masonite  enclosures  met  the  laser  stability  criteria. 

Beam  Diagnostics 

The  distribution  of  laser  beam  intensity  and  fluctuations  in 
intensity  at  the  focal  spot  are  the  two  simplest  observable  qualities 
which  can  be  analytically  predicted  at  this  time.  A large  number  of 
simultaneous  intensity  or  local  power  measurements  are  required  to 
accurately  describe  the  statistics  of  a beam  which  has  been  propagated 
through  a random  medium.  For  a number  of  C.  W.  laser  applications,  the 
time  averaged  laser  beam  power  at  the  focal  spot  is  of  interest  as  well 
as  the  dimension  averaged  spot  size.  A measurement  of  the  time  averaged 
distribution  of  spot  size  will  be  reported  in  this  work. 

To  obtain  the  distribution  of  the  time  averaged  laser  beam 
power,  a circular  aperture  was  traversed  across  the  focal  plane  of  the 
laser  beam.  The  local  beam  power  was  obtained  from  a Spectra  Physics 
Model  401 C power  meter.  The  traverse  of  the  focal  plane  was  accomplished 
by  an  electrically  driven  three-dimensional  traverse  mechanism,  shown  in 
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Figure  11,  located  within  the  light  pipe.  Also  shown  in  Figure  11  is 

-5 

a manual  traverse  which  had  a resolution  of  3.93  x 10  centimeters, 
which  was  discarded  in  favor  of  the  electrical  unit  to  further  reduce 
the  background  noise  level.  The  manual  traverse  exhibited  a higher 
noise  level  due  to  the  required  opening  for  the  manual  traverse  in  the 
light  pipe.  All  traverses  were  made  in  the  vertical  plane  beginning  at 
the  bottom  of  the  focal  spot.  The  focal  plane  was  located  at  a distance 
of  45.7  meters  from  the  test  cabin.  Traverses  were  made  after  initially 
determining  the  location  of  the  focal  plane  along  the  beam  path.  Sub- 
sequent measurements  were  made  after  maximizing  the  power  normal  to  the 

path  of  progagation,  I , I 

max  -^max 

Correlator 

The  scales  of  the  random  fluctuations  were  obtained  from  the  cor- 
relation function  as  will  be  discussed  in  Chapter  IV.  The  correlation 
function  for  temperature  and  velocity  was  obtained  from  a Hewlett 
Packard  Model  3721A  correlator.  The  energy  spectrum  of  the  fluctuations 
was  obtained  by  taking  the  Fourier  transform  of  the  autocorrelation 
function  with  a Hewlett  Packard  Model  3720  spectrum  analyzer.  The  cor- 
relation function  is  given  by: 

limt  ■>  « \ /o  dt  (5) 

where  r is  the  delay  time  between  each  calculated  point  and  A(t)  and 
B(t)  are  the  real  time  functions  for  which  correlation  is  desired.  The 
correlator  calculates  the  correlation  function  for  100  points,  each 
delayed  from  one  another  by  x where  x is  variable  from  1 usee  to  333 
millisec.  The  correlator  accepts  inputs  from  DC  to  250  Kilohertz.  The 
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correlator  is  digital  and  has  two  averaging  processes,  an  arithmetic 
average  or  summation  mode,  and  an  exponential  average.  In  the  summation 
mode  each  point  represents  the  arithmetic  average  of  128  to  128  x 1024 
(N)  calculations. 

In  the  exponential  mode,  each  of  the  100  running  totals  is  calcu- 
lated as: 

Old  Total  + (New  Sample  - Old  Running  Total )/N 
The  number  of  samples  or  sweeps  for  each  point  in  the  exponential  mode 
is  1 to  100  times  N.  A quick  look  at  the  correlation  function  is  pos- 
sible with  the  exponential  mode  of  averaging,  since  a short  time  con- 
stant N x t)  is  used  for  the  early  calculations.  This  allows  the  proper 
amplitude  and  time  delay  increment  t to  be  selected  from  a few  seconds 
of  observation  of  the  correlation  function.  The  probability  density 
and  cross  correlations  were  also  obtained  with  the  same  unit.  The  expo- 
nential averaging  mode  was  used  in  the  experiment  to  be  reported  in  this 
work,  since  it  is  faster  and  suffers  very  little  loss  in  accuracy  to  the 
summation  mode.  The  input  voltage  range  of  the  correlator  is  from  .04 
volts  rms  to  4 volts  rms.  Lower  rms  values  of  temperature  fluctuations 
were  encountered  during  runs  with  zero  heat  addition.  The  TSI  1057 
signal  conditioner  was  used  to  amplify  the  rms  fluctuations  when  below 
the  range  of  the  correlator. 
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CHAPTER  IV 

EXPERIMENTAL  TECHNIQUES  AND  PARAMETER 
CALCULATION  PROCEDURES 


In  order  to  quantitatively  evaluate  the  effect  of  turbulence  on  a 
laser  beam  propagating  through  the  turbulent  medium,  experimental  eval- 
uations of  analytical  expressions  such  as  Sutton's  (Sutton,  1969, 
p.  1741)  equation  for  beam  spread 


n , 1/2 
£ -£r>  i • oL 

a = the  extinction  coefficient 

e = the  angle  of  the  point  where  the  laser  intensity  curve 
has  fallen  to  1/2  Po 


(6) 


a = 2 k2  A <An2>  (1  + i<3i-2.*-2.) 


-5/6 


k = 2tt/X 

An  = the  fluctuations  in  the  index  of  refraction 
A = the  integral  scale  of  turbulence 
Lq  = the  light  path  through  the  turbulent  medium 
D = the  laser  beam  diameter 
a = energy  containing  turbulence  wave  length 
is  typically  required.  The  fluctuating  component  of  the  index  of  refrac- 
tion is  in  turn  related  to  the  fluctuating  component  of  temperature.  In 
order  to  control  and  measure  the  fluctuating  component  of  temperature. 
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both  velocity  and  temperature  must  be  measured.  The  velocity  determines 
the  mean  temperature  rise  of  the  turbulent  medium  through  the  tunnel 
mass  flow  rate.  This  mean  temperature  rise  then  bounds  the  level  in 
fluctuations  in  temperature.  In  the  first  section  of  this  chapter. 

Flow  Field  Measurement,  the  applications,  calibration,  and  frequency 
response  limitations  of  a hot  wire  anemometer  to  measure  both  the  mean 
and  fluctuating  components  of  velocity  and  temperature  will  be  discussed 
in  detail.  The  calculation  of  the  fluctuating  component  of  index  of 
refraction  from  the  temperature  fluctuations  will  be  given  in  the  Flow 
Field  Measurement  section.  Because  the  analytical  calculations  of  laser 
beam  spread  require  the  integration  of  the  fluctuations  in  the  index  of 
refraction  (in  planes  perpendicular  to  the  propagation  direction)  along 
the  laser  propagation  path,  the  temperature,  velocity,  and  fluctuating 
components  of  temperature  and  velocity,  as  well  as  the  calculated  value 
of  index  of  refraction  and  its  fluctuating  component  values,  were 
obtained  along  the  laser  beam  centerlines.  The  procedure  for  calcula- 
tion of  these  mean  profile  values  will  be  given  in  the  Flow  Field  Meas- 
urement section.  The  second  section  of  Chapter  IV,  Correlation  Measure- 
ments, will  discuss  techniques  for  calculating  the  scales  of  turbulence 
which  appear  in  equation  (6).  The  scales  of  turbulence  will  be  obtained 
from  the  autocorrelation  functions  of  temperature  and  velocity.  Section 
two  of  this  chapter  will  formulate  the  autocorrelation  functions  and 
develop  three  calculations  of  the  microscale  and  one  calculation  of  the 
integral  scale.  The  third  section  of  this  chapter  will  discuss  the 
techniques  and  calculations  which  were  employed  to  obtain  measurements 
of  the  laser  beam  spread  in  the  focal  plane,  hj/oQ. 
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Flow  Property  Measurement 

I The  measurement  and  interpretati on  of  velocity  and  temperature  in 

a highly  turbulent  flow  field  is  non-trivial,  even  when  the  flow  in  the 

I / 

test  section  is  very  slow.  The  heat  addition  required  to  create  signif- 
icant fluctuations  in  the  index  of  refraction  results  in  even  larger 
fluctuations  in  both  temperature  and  velocity,  further  complicating  an 
accurate  measurement  of  fluctuating  properties. 

If  Rayleigh  flow  is  considered  through  the  test  section  for  tunnel 

velocities  from  .3  to  18  meters  per  second  and  for  non-dimensional  heat 

-3 

additions  (q/CpT^)  of  3.15  x 10  to  .187,  the  observed  static  and 
stagnation  temperatures  are  found  to  correspond  to  one  another  within 
the  measurement  accuracy  of  temperature.  The  mean  temperature  rise 
ranges  from  .9°  to  53°C.  Static  and  total  pressure  remain  about  (aP/Pq  = 

.0024)  constant  and  equal  to  barometric  pressure.  On  the  other  hand, 
density  and  Mach  number  also  change  significantly  for  Rayleigh  flow  and 
may  not  be  assumed  equal  to  the  preheat  addition  values. 

Even  though  the  use  of  a hot  wire  for  determining  the  flow  and  j 

turbulence  properties  of  the  test  section  gives  only  approximate  results, 
this  method  of  measurement  is  one  of  the  best  and  conventionally  employed 
techniques  available.  The  problems  which  must  be  considered  in  the  use 
of  a hot  wire  for  measurement  of  the  test  section  properties  include: 
separation  of  the  effects  of  temperature,  velocity,  and  density  on  the 
probe,  the  effect  of  high  turbulence  levels  on  heat  transfer,  the  effect 
of  large  values  of  turbulence  on  the  mean  value  of  velocity  and  tempera- 
ture, the  frequency  response  of  the  probe  and  definition  of  the  flow 
regime  (continuum  to  free  molecular)  in  which  the  probe  must  operate. 
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The  flow  regime  is  defined  by  the  Knudsen  number: 


Ku  =1.26/7  M/Re  (7) 

Due  to  the  small  probe  diameter  and  the  wire  operating  temperatures, 
slip  and  even  free  molecular  flow  can  occur  for  very  slow  flows 
(Sandborn,  1972,  p.  55-88).  Knudsen  number  for  the  range  of  conditions 
expected  ranges  from  .27  to  .013.  For  a flow  to  be  considered  to  be 
continuous,  Ku  should  be  £.01.  This  definition  is  somewhat  arbitrary 
and  the  range  of  values  given  above  will  be  considered  to  represent 
continuum  flow.  A change  in  flow  from  continuum  to  slip  flow  should 
exhibit  itself  as  a discontinuity  in  the  heat  transfer  relation. 

The  Stanton  number 

St  = Nu/ (Re  x Pr)  (8) 

is  normally  used  to  describe  fluid  dynamic  heat  transfer;  however, 
(Bradshaw,  1971,  p.  113-114)  there  is  a lack  of  systematic  work  in 
assessing  the  importance  of  the  Prandtl  number  on  wire  calibration,  the 
empirical  relation  given  by  Kramer  and  Van  der  Hegge  Zijnen  gives  a Pr 
dependence  as  follows: 

Nu  = .42  Pr'20  + .57  Pr‘33  Re’50  (9) 

for  a Prandtl  number  of  .714 

Nu  = .393  + .51  Re‘50  (10) 

Bradshaw  indicates  that  Bourke  used 

Nu  = .75  Pr,2°  + .67  Pr'9  Re’50  (11) 

giving: 

Nu  = .70  + .495  Re,3U  (12) 
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Equation  (10)  and  (12)  indicate  that  there  is  a wide  range  of  weak  de- 
pendencies upon  Prandtl  number.  Neither  of  these  efforts  considered 
the  specific  application  to  hot  wires.  The  systematic  investigation  of 
heat  transfer  from  wires  by  Col  1 i s and  Williams  and  the  empirical  rela- 
tionship they  developed  represents  the  most  accurate  correlation  of  hot 
wire  heat  transfer  data  currently  available.  The  heat  transfer  relation 
proposed  (Collis  and  Williams,  1968,  p.  369)  is: 

-.17 


Nu  Op}  ' = A2  + B2  Re11 
Where  Nu  is  the  Nusselt  number  given  by: 

m hd 

Nu  = r 

Where  h,  the  heat  transfer  coefficient,  is  defined  by: 

E2Ri 


5 ■ (17*17)  sh’»'  <T»-Tf> 


(R1+R3 

The  wire  heat  transfer  relationship  then  becomes: 

E2R 


(R1+R3)2  tt  z K (Tw-Tf)  {Tm}  = A2  + B2  Ren 


(13) 


(14) 


(15) 


(16) 


A value  of  .45  is  recommended  for  the  exponent  n (Collis  and  Williams, 
1968,  p.  370)  for  Reynolds  numbers  in  the  range  of  .02  to  44.  In 
this  expression,  Tm  is  a mean  film  temperature  coefficient  defined  as 
(Tw+Tf)/2.  The  thermal  conductivity,  density,  and  viscosity  are  all 
evaluated  at  this  mean  film  temperature,  Tm.  The  wire  temperature  is 
set  by  the  overheat  ratio  selected  for  probe  operation: 

Resistance  Hot 


0 H 


=1+6  (Tw-Tf)  + . 


(17) 


Resistance  Cold 

It  is  clear  that  the  sensitivity  of  the  probe  to  temperature  fluc- 
tuation can  be  increased  by  reducing  the  overheat  ratio  and  likewise  the 
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sensitivity  to  velocity  can  be  increased  by  increasing  the  overheat 
ratio.  One  also  has  the  option  of  operating  a wire  in  a constant  cur- 
rent mode  or  a constant  temperature  mode.  The  constant  temperature  mode 
has  the  advantage  of  a significantly  higher  frequency  response  since  the 
wires  thermal  inertia  is  partially  compensated  for  by  an  electronic 
feedback  loop.  Mean  temperature  and  velocity  can  be  measured  by  opera- 
tion of  two  wires  at  different  overheats  such  that  both  wires  are 
sensitive  to  both  temperature  and  velocity.  This  method  requires  the 
storage  of  large  amounts  of  data  for  future  processing  or  the  instanta- 
neous normalizing,  squaring  at  the  minimum,  and  addition  and  subtraction. 
Some  time  period  must  be  allowed  as  an  averaging  period.  However,  the 
commercially  available  circuits  to  perform  the  arithmetic  functions  were 
found  to  be  unacceptable  for  this  experiment.  The  transfer  functions 
and  frequency  responses  appear  inadequate  for  even  well  behaved  inputs 
such  as  sine  waves.  The  loss  of  experimental  information  due  to  the 
processing  would  have  to  be  defined  with  a known  source  and  this  loss 
then  extrapolated  to  this  experiment.  The  taping  and  later  reduction  of 
data  offers  the  possibility  of  the  best  accuracy;  however,  it  represents 
a significant  undertaking  when  time  samples  less  than  a millisecond  are 
required.  The  experiment  was  expected  to  exhibit  Strouhal  frequencies 
up  to  700  Hertz.  Taped  reduction  of  data  was,  therefore,  considered 
feasible  but  not  attractive  until  an  initial  assessment  of  the  beam 
spread  had  been  made. 

The  approach  which  initially  appeared  to  offer  an  acceptable  meas- 
urement capability  and  accuracy  was  to  operate  a double  parallel  wire, 


one  wire  at  a high  overheat  for  velocity  and  the  other  at  a low  overheat 


for  temperature.  By  operating  the  velocity  wire  in  the  constant  temp- 
erature mode,  and  the  temperature  wire  in  the  constant  mode,  a satis- 
factory selectivity  between  temperature  and  velocity  was  possible. 

A Thermo  Systems  Model  1244  Probe  was  selected  for  measuring  the 
velocity  and  temperature  properties  of  the  tunnel.  The  Model  1244  is 
a double  parallel  3.81  um  diameter  wire  configuration.  The  probe  was 
orientated  with  the  wires  in  the  vertical  plane  in  the  tunnel  to  mini- 
mize convective  and  mutual  interference  effects.  The  wire  was  selected 
over  the  more  rugged  thin  films  due  to  its  increased  frequency  response. 
The  frequency  response  of  the  velocity  wire  of  Probe  1244,  Serial  Num- 
ber 8712,  is  shown  in  Figure  12. 
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Figure  12.  Frequency  Response  of  Probe  8712,  Wire  No.  1, 

Horizontal  Time  Scale  5 x 10'°  Seconds/Division 
Constant  Temperature  - 0 H = 1.5 


The  frequency  response  indicated  by  superposition  of  a 1 KHz  square 
wave  upon  the  bridge  probe  combination,  as  is  recommended  by  TSI  and 
Freymuth  (Freymuth,  1967,  p.  680),  is  given  by: 
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where  t is  the  time  to  return  to  3%  of  the  peak  voltage.  The  frequency 
response  given  by  Figure  12  is  22  KHz  for  the  velocity  wire.  Although 
both  the  film  and  wire  have  adequate  frequency  response  foi  the  velocity 
and  velocity  turbulence  intensity  measurement,  in  the  constant  current 
mode,  which  is  used  to  measure  temperature,  the  electronic  feedback 
circuit  is  removed  and  the  frequency  response  of  the  probe  is  limited 
to  that  naturally  exhibited  by  the  probe  itself.  The  3.81  vim  diameter 
wire,  according  to  the  manufacturer,  has  a 750  Hertz  response  at  the 
-3db  point  in  the  constant  current  mode.  A check  of  the  wire  frequency 
response  may  be  made  by  illuminating  the  constant  current  wire  with  a 
laser  which  is  mechanically  chopped  or  by  using  a pulsed  laser  source. 

A helium  cadmium  laser,  operating  at  5 milliwatts,  was  used  with  two 
mechanical  choppers  to  obtain  the  experimental  frequency  response  of 
the  wire  in  the  constant  current  mode.  A schematic  of  the  arrangement 
which  was  used  to  measure  the  frequency  response  of  the  temperature 
sensor  is  illustrated  in  Figure  13. 


Figure  13.  Schematic  of  Frequency  Response  Measurement 


Figure  14  shows  the  chopped  signal  from  the  Spectra  Physics  401 C power 
meter  located  behind  the  wire,  as  the  upper  trace  with  the  lower  trace 
the  temperature  signal  as  the  probe,  in  still  air,  is  illuminated  by 
the  chopped  laser  beam.  The  power  meter  response  shows  that  the  wire 
is  heated  by  a square  pulse.  The  time  to  decay  to  1/2  VQ  corresponds 
to  812  Hertz,  if  the  rise  time  of  the  pulse  is  neglected.  If  the  time 
to  decay  definition  given  by  equation  (18)  is  used  with  Figure  14,  the 
frequency  response  of  the  temperature  wire  would  be  83  Hertz. 


Figure  14.  Frequency  Response  of  Probe  8712,  Wire  No.  2, 

Horizontal  Time  Scale  ? x lCT-^  Seconds/Division 
Constant  Current  Mode  - 0 H = 1.15 
Chopper  Blade  with  Four  Slots 

Figures  15  and  16  show  a sequence  of  photos  of  the  frequency  response 
of  a wire  in  the  constant  current  mode  obtained  with  an  Infrared  Indus- 
tries Model  830  Chopper,  with  a 40  slot  blade,  as  the  frequency  was 
varied  from  200  to  1200  Hertz.  The  normalized  response  curve  obtained 
from  Figures  15  and  16  is  shown  in  Figure  17.  The  frequency  response 
of  the  single  pulse  of  Figure  14  is  shown  in  Figure  18.  The  frequency 
scale  of  Figure  18  was  expanded  in  Figure  19  to  examine  the  single 
pulse  curve  below  800  Hertz. 
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Figure  17.  Temperature  Frequency  Response  of  Probe  8712  witn 
IR  Industries  Model  830  Chopper  - Multiple  Pulse 
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Figure  18.  Temperature  Frequency  Response  of  Probe  8712  with 
Four-Slotted  Blade  Chopper  - Single  Pulse 


L 


49 


The  analytical  solution  for  the  transient  response  of  a wire  to  a 
step  change  in  fluid  temperature  is  given,  for  first  order  terms,  by 
(Sandborn,  1972,  p.  105): 


4h,d  a 2 -1* 

, f / 4 \ og 4 T 

i‘  P ,C  pCd'  ' wd7^  P C . A1 


K 

- = w 

e w w 


(19) 


w w 


w w 


where  t is  the  time  constant, 
e 

The  first  term  represents  the  contribution  of  the  wire  supports, 
the  second  term  would  be  the  solution  obtained  in  the  limiting  condition 
of  an  infinite  wire  as  1^  ->  0,  and  the  last  term,  the  Joule  heating  of 
8057  inverse  seconds  and  the  support  term  258  inverse  seconds.  The 
approximate  single  pulse,  shown  in  Figure  14,  gives  a 1/t  of  600  inverse 
seconds,  as  the  response  curve  shows  in  Figure  19,  while  the  repetitive 
pulse  shows  a value  of  125  inverse  seconds.  Comparison  of  the  single 
and  multiple  pulse  curves  shows  the  single  pulse  curve  to  be  rather  op- 
timistic *nd  the  wire  is  clearly  not  well  approximated  by  an  infinite 
wi  re. 


It  is  clear  that  any  direct  comparison  of  temperacure  fluctuations 
at  frequencies  in  excess  of  200  Hertz  can  be  subject  to  very  large 
errors.  Although  frequencies  to  approximately  700  Hertz  were  antici- 
pated, the  bulk  of  the  measurements,  which  will  be  reported  for  large 
beam  spread,  exhibit  frequencies  below  100  Hertz. 

In  the  constant  current  mode  the  heat  transfer  relationship  of 
equation  (16)  becomes: 

, -.17  Tm  -.17  El i G 

Nu  {y-1  = f j— 1 jr  s,  K (T  -Tf)  = A2T  + B2T  Re  (20) 

w f 
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Figure  19.  Temperature  Frequency  Response  of  Probe  8712 
for  Single  Pulse 


It  is  desirable  to  select  the  wire  temperature  through  the  overheat 
just  above  the  maximum  temperature  to  be  measured.  For  this  experiment, 
an  overheat,  as  given  by  equation  (17),  of  1.15  was  chosen.  With  the 
wire  in  the  constant  current  mode,  at  an  overheat  of  1.15,  the  experi- 
mental observed  Reynolds  number  dependence  is  shown  in  Figure  20.  Fig- 
ure 20  was  obtained  from  measured  points  with  the  probe  installed  in 
the  tunnel  and  the  least  squares  fitting  of  a straight  line  to  the  data. 
The  least  squares  coefficients  show  that  B^  is  two  orders  of  magnitude 
less  than  Since  the  Reynolds  number  to  the  .45  power  in  this  case 

is  of  order  unity  or  less  the  B^  coefficient  and  Reynolds  number 
dependence  may  be  neglected.  It  is  desirable  for  the  temperature  and 
the  velocity  wires  to  be  located  as  near  one  another  as  possible  and 
still  not  interact  with  one  another.  Mathematically,  the  criteria  for 
neglect  of  convection  was  formulated  (Collis  and  Williams,  1968,  p.  382- 
383)  as: 

Re  > 2 (Gr)1/3  (21) 

where 

g d3  (Tw-Tf) 

Gr  e 2-= r-  (22) 

v Tf 

The  right  hand  side  of  equation  (21)  is  the  order  of  .02  and  probe  wire 
Reynolds  numbers  range  from  .15  to  4.7,  which  justifies  the  neqlect  of 
any  free  convection  interaction  between  the  double  parallel  probes. 

The  radiated  power  is  on  the  order  of  50  microwatts,  which  is  several 
orders  of  magnitude  below  the  probe  power  and  therefore  neqliqible. 
Fiqure  20,  which  varied  the  Reynolds  number  for  the  temperature  wire, 
remains  the  same  regardless  of  whether  the  velocity  wire  is  operating 
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Figure  20.  Reynolds  Number  Dependence  of  the  Temperature 
Wire,  Probe  8712 


at  its  normal  hot  temperature  or  left  at  room  temperature,  providing 
there  is  a measurable  flow  (>1  foot/second)  over  the  probe.  The  free 
convective  and  radiative  effects  are  therefore  neglible. 


The  relationship  obtained  for  heat  transfer  becomes: 

r _ B2T  * 1 K (TW'V  {^1 
IlG 

The  thermal  conductivity  of  the  film  is  given  by: 

K = 2.41  x 10"4  (1  + .00317  Tm  - .0000021  Tm  ) 

where 


(23) 


(24) 


(2b) 


A sample  calibration  curve  for  the  temperature  wire  and  equation  (23) 
is  shown  in  Figure  21. 

By  differentiating  the  relationship  for  temperature  with  respect 
to  time: 


dE 

dt 


(26) 


which  when  divided  by  the  mean  temperature,  the  temperature  turbulence 
intensity  will  be  given  by: 


dT  _ B2j 
T = T " AT" 
2t 


TT  c!E 


IT  + B 


2f 


7T 


E 


(27) 


The  calculation  of  the  index  of  refraction  and  the  fluctuation  in 
the  index  of  refraction,  which  have  been  employed  in  the  majority  of  the 
propagation  efforts,  is  developed  by  Tatarski  as  follows  (Tatarski,  1961 


p.  55-56): 


n = 1 + 79  x 


(20) 


where  P„  is  pressure  in  millibars  and  T is  the  temperature  in  °K,  provid- 

D 

ing  the  effect  of  humidity  on  the  value  of  index  of  refraction  is 
neglected.  Differentiation  of  the  index  of  refraction  with  respect  to 
time  yields : 


= 79 
dt 


x 10 


-6 


f _ 
vdtT 


dT  i 
dtT^ 


(29) 


with  the  assumption  of  a locally  isentropic  process 

dT  _ y-1  dp 
T “ y P 

at = 79  x 10~6  {m x p " 1)} 

with 


1.4 


Q = 1.975 
dt 


x 10 


-4  PdT 
dtT2 


(30) 

(31) 


(32) 


After  the  temperature  of  the  fluid  has  been  determined  by  the  wire 
operating  in  the  constant  current  mode,  the  measurement  of  velocity  may 
be  undertaken.  The  heat  transfer  relation  with  the  wire  operating  in 
the  constant  temperature  mode  was  given  by  equation  (20).  This  relation 
was  programmed  to  least  squares  fit  a curve  to  the  calibration  data  and 
calculate  the  curves  least  square  coefficients.  The  zero  velocity 
point  has  been  neglected  in  fitting  the  curve  to  the  calibration  data. 
The  zero  velocity  point  has  been  neglected  since  there  is  only  free 
convection  heat  transfer  with  neglible  radiation  at  zero  velocity 
(Bradshaw,  1971,  p.  116).  Calibration  curves  were  also  plotted  using 
the  zero  flow  point;  however,  this  results  in  larger  errors  for  all 
calibration  points.  Sample  calibration  points  are  shown  in  Figure  22. 
Figure  23  shows  a least  squares  fit  of  calibration  data.  The  calibration 
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curve  shows  no  free  convection  effects  at  low  values  of  Reynolds  number 
and  shows  no  transition  from  slip  flow  to  continuum  flow  over  the  cali- 
bration range.  The  fluctuations  in  velocity  may  now  be  calculated  and 
the  turbulence  intensity,  Tu,  obtained: 

Tu  = (33) 

^ {Nu  {^1  17}  = ^ (A2  + B2  Re’4b)  (34) 

Rewriting  to  collect  the  temperature  effects  in  one  term: 

5t  (E?)  = ft  11  1 K (Tw-Tf)  ^“}17}  (A+B  Re‘45)  (35) 

Let  the  first  bracket  be  defined  as  6 (T^). 

Giving  the  density  used  in  Reynolds  number  calculations  as: 

ps  PG  Ts  mm 

" TyrfTT7  (3fc) 

2 

Viscosity  in  the  Reynolds  number  is  likewise  evaluated  at  the  mean  film 


temperature : 


T,  ,+Tr  3/2 

/ W i \ / ioo 

M " uo  (~273~'  (Tw+V 


+ 113 


with  the  temperature  in  °K. 

Differentiating  the  expression  for  the  square  of  the  wire  voltage  with 
respect  to  time  yields: 

2E  ^ ' .45  B2T  iT’55  B(Tf)  + ^ (B2T  Re’45  B(Tf))  + 


dt  ^A2T  6^Tf^ 
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While  evaluation  of  the  temperature  contributions  could  be  carried  out, 

there  is  a very  complex  relationship  for  an  accurate  description.  As- 

Tm  * 

suming  viscosity,  thermal  conductivity,  and  { j — } constant  would 

result  in  the  following  expression  for  velocity  turbulence  intensity: 

T du  _ 2E_dE  _ 
u ~ U .45  Re-^  6(Tf) 

r TW  . J_  . A2T  , 

1 (Tw+Tf ) • 55  (Tw+Tf)  * .45  B2T  Re-4:>  .45  dTf  (39) 

Tw  - Tf 


The  temperature  term  corrections,  as  given  by  equation  (39)  were  included 
for  calculation  of  Tu  values  on  the  tunnel  centerline.  For  the  case  of 
zero  heat  addition,  the  term: 


Tu 


? e dE 

. 45  Re • u b B2T  K h a“(R  (Tw -T f ) 


(40) 


was  calculated  for  the  traverses.  This  relationship  is  shown  in  Figure 
24  for  the  wire  calibration  curve  of  Figure  23. 

Because  the  levels  of  turbulence  are  exceeding  high  and  beyond  the 
range  of  values  in  the  conventional  appl  ication^f  wires,  the  effect  of 
high  levels  should  be  examined  on  both  the  heat  transfer  and  also  the 
possible  non-linear  effects  on  determination  of  mean  values.  Sandborn 
(Sandborn,  1971,  p.  235,  Fig.  6.6)  shows  significant  increases  in  heat 
transfer,  1.65  times,  for  a Tu  of  12%  at  A/d  ratios  of  1.5.  This  effect 
is  reduced  to  the  order  1.1  times  the  heat  transfer  at  a A/d  of  8,  where 
• is  the  scale  of  turbulence  and  d the  hot  wire  diameter.  Since  the 
ratio  of  /d  for  this  experiment  is  of  the  order  120,  this  effect  should 
• negligible  High  levels  of  turbulence  can  also  cause  the  mean  value 
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Figure  24.  Tu  for  Probe  8712  at  a Fluid  Temperature  of  23°C 


of  voltage  observed  on  the  probe  not  to  correspond  to  the  mean  value  of 
velocity.  When  this  happens,  it  is  necessary  to  use  more  elaborate 
techniques,  either  analytic  or  graphical,  to  make  corrections  to  the 
mean  value  of  velocity.  For  evaluation  purposes,  the  graphical  tech- 

[ ! 

nique  was  used  to  determine  the  highest  turbulence  intensities  observed 
and  their  effect  on  the  determination  of  the  mean  velocity.  This  prob- 
lem arises  when  the  non-linear  region  of  the  calibration  curve  is  used 
as  illustrated  in  Figure  25.  If  the  rms  value  of  the  observed  fluctua- 
tion is  superimposed  symmetrically  upon  the  measured  mean  value  of  D.C. 
voltage  as  shown:  (V-j  + may  not  represent  the  mean  value  of  veloc- 

ity implied  by  the  measured  mean  value  of  voltage,  for  turbulence 
intensities  of  the  order  of  20%.  There  was  no  observed  difference  in 
mean  value  of  velocity  at  the  lowest  velocity,  where  this  consideration 

I 

is  most  critical.  For  a 60%  turbulence  intensity  level,  the  correction 
amounted  to  a 7%  error  in  Vm  at  5 feet  per  second  and  a subsequent  3% 
error  in  Tu.  A non-linear  correction  to  Vm  was  therefore  not  considered 
to  be  required  in  data  reduction. 


VI  V2 
Vrn 


VELOCITY 

Figure  25.  Schematic  of  Graphical  Computation 
of  Vm  for  Large  Values  of  Tu 
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Tunnel  profiles  of  mean  velo<  ity,  i . 

and  rms  temperature  traverses  cn  the  horizontal  centerline  of  : i , 

were  made  at  tne  three  i_/Dp  stations  for  which  v.e  beam  spread 
tained.  A traverse  was  also  performed  in  t!.-  vcati  • ; -Ur. 
that  any  natural  convective  effects  were  small  The  hea  er  war  , 
we  re  045  centimeter  diameter  Ni chrome  wires,  for  which  c - .1 

oer  s of  the  order  of  2.37.  A Reynolds  number  of  ,'ess  tr.v 
be  required  to  neglect  free  convective  effects.  Free  cc/nve." 
may  be  neglected  since  the  Reynolds  number  for  tne  tunri  i ,•  ; ic- 

on the  order  of  19  to  560,  which  satisfies  the  cm  -1 

For  each  incremental  Nusse’.t  number  and  Reynol-  ;e 
was  obtained  which  consisted  of  the  mean  velocity,  temperature  r . . 
locity,  end  rms  temperature.  The  nominal  sen  ence  <,“ii 
beam  traverses  was  no  heat  addition,  low  heat  (527  watt.;, 

(783  watts),  high  heat  (1166  watts),  arid  then  a rerun  at  no  he  c i 
tion  to  check  that  the  tunnel  mean  velocity  had  remained  ronstaut.  Tit: 
traverse  extends  to  within  2.9  centimeters  of  the  tunnel  side  wells, 
two  wire  probe  holders  prevented  closer  approach  of  th 

Traverses  were  performed  at  four  L/Dp  ratios,  three  of  > ; ' ,h  are 
illustrated  in  Figure  26.  For  each  traverse,  three  core.1  rt  f vs 

locity  and  temperature  were  obtained  at  i r to r 01  ■ -c  A 

and  C, . These  stations  are  symmetric  about  he  to: 
are  shewn  schematically  in  Figure  27  Twc  'lewlet  Pa  d 
i^.ters  were  used  for  me ssurament  of  the  flue  uatioi  mp  r r v 
ity  and  temperature.  The  rms  meter  plotter  combination  was  ^ .1 

witn  3 00  Hertz  sine  wave 

6 i 
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Tunnel  Profile  Calculations  for  Mean  Temperature,  Mean  Velocity, 


Mean  Index  of  Refraction,  RMS  Temperature,  RMS  Velocity,  RMS  Index  of 
Refraction,  Tu,  Tj,  and  Tn 

For  each  profile  the  centerline  value  of  temperature  is  first  cal- 
culated, then  a mean  value  of  temperature  is  determined  by  continuously 
digitizing  the  temperature  with  T(Ti  + Ti+-|)/2]  ax^  accumulated.  In 
excess  of  a 100  data  points  is  accumulated  in  the  calculation  of  a mean 
value.  In  most  cases  the  centerline  value  is  representative  of  the 


value  determined  by: 


T.  + T. . , 


))  AX1 


Knowing  a mean  value  of  temperature  next  allows  a mean  value  of  velocity 
to  be  calculated  on  the  tunnel  centerline.  A centerline  velocity  is 
calculated  followed  by  the  determination  of  a mean  value  of  velocity  by 
a continuous  digitization  of  the  velocity  traverse  similar  to  the  temp- 
erature calculation,  doing  a simple  trapezoid  integration,  and  then 
dividing  the  result  by  the  accumulated  machine  traverse: 

U.  + IK. 

Z(— 5 ^—4  Axi 

^ / a n \ 


A minimum  of  75  points  would  typically  be  accumulated  in  calculation  of 
U.  The  centerline  and  mean  value  again  generally  agree.  A mean  Reynolds 
number  is  also  calculated. 

The  fluctuation  in  temperature  is  next  calculated: 

dT.  + dT.., 

z(— — 5 -)  Axi 

AT  - / A 1 \ 
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and  combined  with  the  mean  temperature  to  give  the  temperature  turbu- 
lence intensity.  The  index  of  refraction  and  the  fluctuation  in  the 
index  of  refraction  is  next  calculated,  assuming  that  the  pressure  is 
barometric  pressure  as  was  indicated  by  Rayleigh  flow  considerations. 
The  integral  of  the  index  of  refraction  and  the  square  of  the  fluctu- 
ating components  are  to  be  later  used  in  calculation  of  beam  spread,  so 
these  values  were  also  printed  out.  The  rms  level  of  both  temperature 
and  velocity  show  considerable  variation  and  so  a maximum  and  minimum 
value  of  temperature  turbulence  intensity  and  velocity  Tu  is  also 
calculated.  Approximately  300  points  are  typically  accumulated  in  the 
measurement  of  a mean  fluctuating  value. 

Finally,  the  mean  value  of  velocity  turbulence  intensity  is  ob- 
tained in  a similar  fashion  after  calculating  a mean  value  of  the  rms 
level  of  velocity: 


dU  = 


dll.  + dU . , , 
z( — ^—5 -+-)  Axi 


2 

Z Axi 


(44) 


Examination  of  the  profiles  shows  that  both  temperature  and  veloc- 
ity exhibit  a turbulent  distribution.  Both  the  mean  temperature  and 
velocity  are  reasonably  uniform  across  the  traverse,  with  the  center- 
line  values  being  representative  of  the  spatial  mean  values.  When  beam 
spread  measurements  were  made,  the  centerline  tunnel  velocity  was  set 
to  the  same  centerline  profile  conditions. 

The  rms  profiles  showed  much  larger  spatial  variations  titan  were 
anticipated  and,  therefore,  a space  average  value  is  very  important 
for  the  selection  of  a realistic  rms  average  value  and  turbulence 
intensities. 
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Determination  of  Integral  Scale  and  Microscale  of  Turbulence  from 
Correl a t i on  Functions 

The  prediction  of  beam  spread  by  optical  calculation  requires  in- 
formation about  the  integral  scale  and  microscale  of  turbulence  for  the 
flow  field  through  which  the  beam  is  propagating.  These  scales  of 
turbulence  may  be  deduced  from  the  correlation  function  of  velocity, 
temperature,  and  optical  fluctuations.  A Hewlett  Packard  3 7 2 1 A corre- 
lator was  employed  to  obtain  the  autocorrelation  function  of  interest 
in  the  time  domain.  The  algorithm  of  the  correlator's  computer  was 
periodically  verified  by  inputing  a sine  wave  of  known  amplitude  and 
frequency. 

In  this  section  the  correlation  function  will  be  defined  and 
Taylor's  Hypothesis  will  be  invoked  in  order  to  transform  the  results 
from  the  time  domain  to  the  spatial  domain.  The  integral  scale  of 
turbulence  will  then  be  obtained  in  a straight  forward  manner  from  the 
correlation  function.  The  microscale  of  turbulence  may  be  deduced  from 
the  correlation  function  by  several  procedures.  These  procedures  will 
be  discussed  and  the  microscale  will  then  be  determined. 

The  integral  scale  of  turbulence  represents  the  larger  set  of 
eddies  contained  in  the  flow  field.  These  eddies  are  controlled  by  the 
boundary  conditions  and  are  proportional  to  the  greatest  distance  be- 
tween two  points  in  the  flow  field  for  which  there  exists  a correlation. 

The  correlation  function  may  be  obtained  for  fluctuation  in  both 
scalar  and  vector  fields.  The  correlation  function  for  temperature 
(scalar  field)  and  velocity  (vector  field)  are  shown  schematically  in 
Figure  28  with  the  mean  flow  in  the  X direction. 
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Figure  28.  Schematic  of  Longitudinal  and  Lateral 
Correlation  Functions  (Hinze,  1959, 
p.  30,  Fig.  1-14) 

Here  f(r)  and  g(r)  are  respectively  the  longitudinal  and  the  lateral 
velocity  correlation  functions  and  T(r)  is  the  temperature  correlation 
function.  The  lateral  correlation  function  will  be  the  experimental 
quality  to  be  measured  to  determine  the  scales  of  turbulence.  The 
correlation  function  is  a fundamental  property  which  is  used  to  charac- 
terize random  processes.  If  M is  denoted  by  the  distance  c and  N by 
C + x^,  then  the  lateral  correlation  function  is  given  by: 

9(r)  = U2  U)  U2  U + x2)  ^ (45) 

g(r)  = [U2  (c)  - U2c]  [U2  (c  + x2)  - U2  U + x2)]  (46) 

The  integral  scale  of  turbulence  can  now  be  defined  in  terms  of  the 
correlation  function  by: 

A = r g(r)  dx  (4/) 

o 

Following  Tatarski  (Tatarski,  1961,  p.  29),  the  qualitative  connection 
between  the  microscale  and  the  integral  scale,  which  was  stated  in 
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equation  (2),  will  be  outlined.  The  turbulent  kinetic  energy  (u'-u1), 
which  is  converted  to  thermal  energy  by  dissipation,  can  be  expressed 
by: 

z ^(u1 -u' )/x  (4bj 

where  t is  a characteristic  period  defined  as  t = for  a fluctuation 
or 

e ^(u’-u'J/t  = u 1 3/A  (49) 

This  is  assumed  to  hold  for  all  eddies  except  the  smallest  for  which 

e 'v  v 1)  2/£  2 (50) 

0 0 

and 

u.  - (•  <„>,/3  i5’) 


or 


&0  ^ rV3/e 


(52) 


where  zQ  is  the  microscale  of  turbulence.  The  microscale  can  then  be 
expressed  in  terms  of  the  integral  scale  by 
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(ReA) 


3/4 


(53) 


At  this  point  it  should  be  noted  that  the  prediction  of  the  scales  of 
turbulence  and  relationship  between  the  scales  is  in  reality  very 
empirical.  Previously  equation  (2)  was  introduced  which  gave  a some- 
what different  empirical  relationship  between  £o  and  A than  equation  (53). 
Corrsin  gives  equation  (2)  as: 
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where 


Re. 


Re 


Re . ^ /Re 

£o 


(55) 

(56) 


assuming  A/Dp  and  u ' /II  are  independent  of  Reynolds  number  at  a fixed 
distance  downstream  of  a turbulence  generator  such  as  a grid.  Then  A^ 
is  assumed  to  be  constant.  The  value  of  A ranges  from  2 to  30  for 
Biplanar  cylindrical  grids  (Corrsin,  1963,  p.  528)  (Goldstein,  1965, 
p.  228).  Although  the  deduction  of  a microscale  and  an  integral  scale 
of  turbulence  from  the  autocorrelation  function  is  recognized  as  an 
inadequate  specification  of  actual  scale  sizes  for  many  turbulent  flow 
fields,  these  representative  scales  are  employed  in  many  optical  cal- 
culations such  as  those  of  equation  (6)  for  beam  spread. 

There  exist  a variety  of  ways  to  deduce  the  microscale  of  turbu- 
lence from  the  correlation  function,  three  of  which  will  follow.  Be- 
cause the  flow  is  considered  to  be  homogeneous,  g(x^)  = gt^),  Hinze 
(Hinze,  1959,  p.  35-36)  shows  that 
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for  n + m odd  and  n + m even  by: 
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(58) 


x^O 


The  general  and  complete  shape  of  the  correlation  function  of  course 
depends  upon  the  nature  of  the  turbulent  flow  field.  However,  at  x^-O, 
the  correlation  function  can  be  described  in  terms  of  its  derivatives 
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at  the  origin.  Expanding  g ( x ) in  a Taylor  series,  the  noting  that  g(x) 
is  symmetrical  with  respect  to  x.  : 


gU, 


+ 1_  (x  \2  ri22  i 


3x,/Jx-  = 0 


+ . 


(59) 


Sjbst’tut;on  of  the  partial  derivative  of  g into  the  expansion  gives: 

(x?)2  1 3U-.  2 

»'*?>  ’ 1 - n (in*  . ♦ • • • (»> 

• u ” u 

which  for  srrial'1  values  of  may  be  approximated  by  a parabolic  function 


of  x^: 


g(*2)  * 1 


(61) 


or 

I ] | y i . J__  r 3U1  -i2 

fc,  -v  2 L?xfJx0  = 0 2TuTT  ^x,-1 

u 2 2 1 2 


(62) 


r-  which  ? is  re  microscale  of  turbulence, 
o 

The  above  expressions  for  the  microscale  suggest  two  techniques 
for  ct  lasurement  of  the  microscale,  g ( ) near  the  origin  can  be 
■s-ed  to  ; t a parabola  through  tne  origin  and  another  point  near  the 
origin. 
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Figure  29.  Parabolic  Definition  of  the 
•Microscale  of  Turbulence 
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The  second  expression  of  equation  (62)  can  be  evaluated  either  with 
electronic  differentiation  of  the  rms  value  squared  of  the  parameter 
from  which  the  correlation  function  is  calculated  or  by  mathematically 
taking  the  second  partial  of  the  autocorrelation  curve  at  = U.  The 
frequency  response  and  average  time  for  the  measurements,  as  well  as 
the  response  to  a random  wave  function  are  questionable  using  electronic 
differentiation.  The  electronic  differentiation  was  discarded  in  favor 
of  fitting  a cubic  equation 

g(x2)  = Mq  + M-|  x + M2x2  + M3x3  (63) 

to  the  first  four  points  of  the  autocorrelation  function  and  then 
solving  for  a2  g(x2)/9x22. 

A third  method  of  calculation  of  the  microscale  may  be  developed 
by  examination  of  the  predicted  behavior  of  the  autocorrelation  curve 
as  x2  becomes  large.  To  accomplish  this  calculation  procedure,  it  is 
necessary  to  relate  the  correlation  function  to  the  structure  function. 
For  the  measurements  to  be  reported,  the  turbulent  flow  field  may  be 
considered  to  be  locally  homogeneous  (Batchelor,  1953,  p.  133)  at  any 
given  transverse  station.  The  assumption  of  homogeneity  of  the  turbu- 
lent flow  field  was  confirmed  with  x probe  hot  wire  measurements  for 
both  velocity  and  temperature  at  an  L/Dp  tunnel  station  of  28. 

A homogeneous  flow  field  at  a given  tunnel  station  may  be  repre- 
sented by  stationary  process.  Following  Tatarski  (Tatarski,  1961, 
p.  8-33),  the  correlation  function  is  represented  in  the  time  domain 
(all  scale  calculations  will  subsequently  be  converted  to  the  time 
domain)  approximately  by  a random  function  with  stationary  increments 


72 


by  examining  the  function  about  some  small  increment  ir  time,  t.  The 


stationary  function  is  obtained  by  subtracting  the  random  function 


F ( t ) . 


Ft  (x)  = f u + x)  - f (0 


Then  employing  the  identity 


(a-b)  (c-d)  = \ [(a-d)2  + (b-c)2  - (a-c)2  - (b-d)2]  (65) 

The  stationary  correlation  function  becomes 

g Ur?2)  = g U-|)  g U2)  (66) 

g (t-|,t2)  = \ [g( c-j+x)  - g(t2)]?-  + ^ [g( ) - gU2+x)]2  (67) 

g (trt2)  = - \ [g(t-|+x)  - g(c2+x)]2  - ^ [g( c-, ) - g(c2)]:  (68) 

This  linear  combination  of  stationary  functions  is  defined  as  the 


structure  function 


Dg(x) 


or  for  a random  process  with  a zero  mean  value 


y«)  - 1 

Since  g(t)  is  stationary 


i-xT]2  + g(t)]2  - 2gU+x)  g(c) 


which  gives 


[g(t)]2  = [g(c+x)]?  = g ( 0 ) 


Dg(x)  = 2 [g(0)  - g(x)] 


There  is  a transverse  and  longitudinal  structure  function  which  are 
functionally  related. 

For  x such  that  n <<x  <<A  dimensional  arguments  give  the  "two/ 
thirds  law"  for  the  form  of  the  structure  function: 


Dg(x)  = c^x)2'3  (73) 

for  values  of  x„  less  than  a 
2 o 

Dg(x)  = c2(x)2  (74) 

Although  the  constants  vary,  the  form  and  functional  dependence  for  both 
the  longitudinal  and  lateral  structure  function  is  the  same  as  given  in 
equations  (73)  and  (74).  Recalling  that  the  structure  function  was 
given  in  terms  of  the  correlation  function  by  equation  (72),  the  struc- 
ture function  may  be  obtained  by  displacing  the  origin  to  the  first 
point  of  the  autocorrelation  function,  reflecting  the  curve  about  the 
new  axis,  and  multiplying  the  ordinate  by  2.  The  microscale  may  be 
found  by  equating  So  - S = 0 or  x = So/S.  c,  of  equation  (73)  may 

X X I 

be  calculated  at  a time  greater  than  the  time  of  the  point  of  inflection 
for  the  autocorrelation  curve.  The  microscale  may  then  be  taken  as  the 
point  where  the  cubic  is  less  than  the  value  of  the  autocorrelation 


curve  as  is  shown  in  Figure  30. 


C*  X 


Figure  30.  The  Scales  of  Turbulence  from  the  "Two/Thirds  Law" 
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The  integral  scale.  A,  is  the  value  at  which  the  curve  for  the  structure 
function  begins  to  saturate  as  is  shown  in  Figure  30. 

All  turbulence  scale  measurement  techniques  are  more  convenient  to 
employ  if  the  averaging  process  is  accomplished  in  the  time  domain 
rather  than  the  space  domain  for  which  the  correlation  functions  were 
originally  defined.  If  the  time  varying  correlation  function  is  exam- 
ined at  a single  point  in  the  flow  field,  similar  relations  are  obtained 
for  the  expansions  of  the  correlation  function  (Hinze,  1959,  p.  40-41). 

9(t)  = 1 - £ (75) 

l l 9ul  ° 

^=2^[3F-]2  <76> 

To  ^U1  9t  t=0 


with  an  integral  time  scale  given  by: 

tq  = dt  g(t)  (78) 

The  exact  relationship  between  g(t),  f(x),  and  g(x)  is  unknown  (Hinze, 
1959,  p.  40-41).  If  the  turbulent  field  is  homogeneous  with  a constant 
mean  velocity  in  the  x-j  direction  = IF,  then  the  temporal  and  spacial 
observations  may  be  interchanged  through  Taylor's  hypothesis: 

ft  = " ^ fx  (79> 

providing  u-j  is  small  compared  to  U.  The  requirement  for  u^/lF  results 
from  consideration  of  the  equation  of  motion: 

t ,ru,  , 8ui  . 1 sip  . v8;“) 

dt  dt  ^ ui  p 3x^  3x|(3xk 
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Since  p-j  « pu^  the  approximation,  "Taylor's  Hypothesis,"  will  be  valid 
if  u-j/U  <<1.  Taylor's  Hypothesis  will  be  applied  in  this  case  also  when 
the  u^/U  criteria  is  not  met.  The  connection  between  the  integral  time 
scale,  t,  and  the  integral  spatial  scale  is  given  by: 

Ajj  = U Ty  = U f g(ty)  dt  (81) 
At  = U tt  = U /“  g(tT)  dt  (82) 
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Sample  autocorrelation  functions  will  next  be  shown  to  illustrate 
the  range  of  conditions  encountered  in  this  experiment.  Approximately 
1700  autocorrelations  of  temperature  and  velocity  were  taken  to  document 
the  turbulent  flow  field.  Three  sample  curves  are  shown  in  Figures  31, 
32,  and  33.  The  upper  most  curve  is  the  autocorrelation  function  and 
the  lower  curve  the  Fourier  transform  of  the  autocorrelation  function 
which  is  the  energy/frequency  spectrum  of  the  turbulence.  Figures  31 
and  33  are  velocity  autocorrelations  and  Figure  32  is  a temperature 
autocorrelation  at  the  centerline  of  the  tunnel  at  the  L/ Dp  station  of 
45.  Comparison  of  Figures  31,  32,  and  33  with  Townsend's  (Townsend, 
1956,  p.  4-20)  examples  would  indicate  that  Figure  31  is  representa- 
tive of  isotropic  turbulence  with  a single  eddy  size.  Figure  32  appears 
to  be  the  sum  of  two  distributions  which  indicates  a field  of  isotropic 
turbulence  with  two  characteristic  eddy  sizes,  and  Figure  33  has  a 
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Figure  31.  Autocorrelation  of  Velocity  - Single 
Scale  of  Turbulence 


Figure  32.  Autocorrelation  of  Temperature  - Two 
Scales  of  Turbulence 


Figure  33.  Autocorrelation  of  Velocity  - Multiple 
Scales  of  Turbulence 
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steep  decay  near  origin  indicating  isotropic  turbulence  with  eddies  of 

many  differing  sizes.  The  frequency  spectrum  in  each  use  indicates  a 

similar  conclusion.  The  multiple  scales  can  cause  difficulties  in  curve 
2 

fitting  for  the  method  of  microscale  measurement  and  near  the  origin 
there  are  clearly  numerous  points  which  have  a quadratic  function  de- 
pendence with  the  abscissa,  time. 

Laser  Beam  Power  Measurement  and  Calculation 

The  spread  of  the  beam  was  measured  as  an  intensity  profile  in  the 
focal  plane.  This  section  will  discuss  the  linearity  and  frequency 
response  of  the  optical  detector,  the  selection  of  an  aperture  for  the 
focal  plane  traverse,  and  finally  the  calculation  of  the  laser  power 
profile  and  wander  of  the  beam  centroid  from  the  measured  traverse. 

Detector  Calibration  and  Frequency  Response 

The  detector  used  behind  the  aperture  was  a Spectra  Physics 
Model  401C  power  meter.  The  power  meter  was  found  to  have  a linear 
response,  +2%.  The  frequency  response  of  the  power  meter  was  measured 
with  rotating  mechanical  choppers  and  the  helium  neon  laser.  The  fre- 
quency response  of  the  power  meter  is  shown  in  Figure  34  where  the  power 
meter  frequency  was  taken  as  the  average  of  four  to  nine  cycles.  Figure 
35  shows  that  frequency  is  faithfully  tracked  to  25  Kilohertz.  The 
meter  output  was  foun<t  to  be  constant  to  3000  Hertz  with  a drop  to  52% 
of  the  zero  frequency  output  at  20  Kilohertz  as  shown  in  Figure  35.  The 
corresponding  frequencies  associated  with  the  integral  and  microscales 
of  turbulence  in  this  experiment  were  generally  below  100  Hertz.  The 
Fourier  transform  of  the  autocorrelation  functions  of  velocity  and 
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PCWER  METER  FREOJENCY  IN  HERTZ 


Figure  34.  Power  Meter  Frequency  Response 
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temperature  also  confirmed  the  fact  that  there  is  very  little  turbulent 
energy  at  the  higher  frequencies.  The  frequency  response  of  the  detector 
was,  therefore,  considered  adequate  for  this  experiment. 

The  Selection  of  the  Scanning  Aperture  Size  at  the  Focal  Spot 

A number  of  methods  may  be  employed  to  examine  the  intensity 
distribution  of  a laser  beam  in  the  focal  plane.  Microdensitometers, 
videcons,  and  apertures  have  been  successfully  used.  Apertures  are 
simple  and  inexpensive  if  the  integration  of  the  non-uniform  distribu- 
tion can  be  performed  and  the  traversed  distribution  be  shown  to  be  an 
accurate  representation.  It  would  appear  that  the  smallest  aperture 
would  give  the  most  accurate  measurement  of  the  intensity  distribution. 
However,  the  smaller  the  aperture,  the  more  sensitive  it  is  to  scintil- 
lation effects.  The  signal  is  also  reduced  as  the  aperture  diameter  is 
decreased,  while  the  background  noise  level,  which  includes  building  and 
room  vibrations  and  circulation  in  the  external  optical  path  remains  at 
a constant  level,  independent  of  aperture  size.  As  a result,  as  smaller 
apertures  are  used,  the  signal  to  noise  level  is  proportionally  reduced. 
One  hundred  (100)  to  800  micron  diameter  apertures  were  evaluated  to 
choose  an  aperture  diameter  for  the  traverses.  The  aperture  percent  of 
beam  diameter  in  the  focal  plane  has  been  tabulated  in  Table  1 for  the 
lasers  and  beam  diameters  to  be  reported  in  this  work.  Table  1 shows 
that  the  aperture  diameter  percent  of  beam  diameter  ranged  from  5.5%  to 
47%.  Figure  36  shows  a comparison  of  the  results  of  profiles  obtained 
with  200  micron,  400  micron,  and  600  micron  apertures  compared  to  a 
gaussian  distribution.  The  400  micron  aperture  with  the  50  millimeter 
diameter  helium  cadmium  beam  is  observed  from  Figure  34  to  track  the 


81 


°I  / I 


82 


Figure  36.  Comparison  of  Focal  Plane  Aperture  Intensi 
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gaussian  distribution  to  within  2%.  Table  1 shows  that  the  400  micron 
aperture  with  the  50  millimeter  diameter  helium  cadmium  beam  has  an 
aperture  to  focal  spot  ratio  of  23%. 

The  other  aperture  diameters  evaluated  exhibit  comparatively 
poorer  tracking  ability  as  shown  in  Figure  36.  The  ratio  of  aperture 
diameter  to  focal  spot  diameter  should  be  maintained  constant  in  com- 
paring helium  neon  and  helium  cadmium  laser  beam  spread  results.  For 
example,  a 600  micron  diameter  aperture  would  be  selected  for  a helium 
neon  beam  of  50  millimeters  diameter  and  a 400  micron  aperture  would  be 
selected  for  a 50  millimeter  diameter  helium  cadmium  beam. 

The  Calculation  of  Integral  Laser  Power 

The  power  profile  obtained  behind  the  circular  aperture  was 
digitized  and  the  profile  integrated  on  a Hewlett  Packard  9800  series 
digitizer  and  calculator.  The  focal  spot  is  located  prior  to  taking  a 
set  of  data.  Prior  to  running  a profile,  the  intensity  or  power  is 
maximized  in  the  plane  of  the  focal  spot.  A single  traverse  requires 
30-60  seconds  and,  therefore,  includes  a time  average  of  the  focal  spot 
motion. 

The  integration  calculation  assumes  that  the  profile  is  sym- 
metric. The  original  experimental  traverses  were  obtained  on  27.9  x 
43.2  centimeter  paper.  Due  to  the  large  size  of  the  trace,  no  correc- 
tions were  made  for  rotation  of  the  trace  relative  to  the  digitizer 
table  during  the  integral  calculation.  A base  line  was  recorded  on  each 
trace  and  a line  parallel  to  the  base  line  was  aligned  with  the  digitizer 
base  line.  Vertical  and  horizontal  scales  are  entered  by  digitizing  the 
peak  intensity  and  the  two  end  points  of  the  laser  profile.  Because  the 
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Table  1.  Aperture  Percentage  of  Beam  Diameter 
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Figure  38.  14.03  mm  Diameter  He  Cd  Laser  Beam  Traverse  at  an 

L/Dp  Station  of  28,  Zero  Heat 
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centerline  of  the  focal  spot  is  moving,  the  average  centerline  is  ini- 
tially estimated  and  the  zero  for  the  digitizer  shifted  to  this  value. 
Successive  pairs  of  points  at  the  same  intensity  are  entered  and  a mean 
radius  and  centerline  location  are  calculated.  Two  sample  intensity 
profiles  are  shown  in  Figures  37  and  38.  Figure  38  is  the  zero  heat 
case  and  Figure  37  represents  a medium  heat  case  for  a 14.03  millimeter 
helium  cadmium  beam  at  an  L/Dp  station  of  28.  The  intensity  or  power 
profile  is  reduced  to  a single  average  half-width  power  profile.  The 
area  under  the  fit  of  a straight  line  through  pairs  of  successive  points 
is  integrated  and  then  revolved  through  2tt.  The  integral  under  the 
straight  line  approximation  is  accumulated  as  total  power  as  a function 
of  the  radial  position,  assuming  the  average  half-width  distributions  to 
be  symmetric  in  the  focal  plane.  The  actual  values  of  centerline,  aver- 
age radius  and  incremental  accumulated  power  are  pointed  out.  For  plot- 
ting the  results,  the  power  is  normalized  to  the  no  flow  value  and  the 
average  half-width  to  the  value  of  radius  at  the  no  flow  integrated 
power  value  point.  The  digitized  integral  curves  of  Figures  37  and  38 
are  shown  in  Figures  39  and  40. 

The  volume  of  a right  circular  cylinder  was  calculated  by  hand 
and  with  the  program  to  check  the  calculation  procedure  and  assure  that 
rotation  of  the  records  was  not  significant.  The  hand  calculation 
agreed  within  .017%  of  the  machine  calculation  for  eight  increments. 

The  zero  flow,  zero  heat  calculations  typically  involve  approximately 
60  pairs  of  points.  The  digitizer  has  a machine  resolution  of  .0254 
machine  centimeters.  The  centerline  is  shifted  to  the  peak  value  of  the 
laser  intensity  distribution  at  the  beginning  of  the  integral  laser 
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Figure  40.  Integral  Laser  Power,  Medium  Heat 
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power  calculation.  The  shift  of  the  centerline  location  is  convenient 
for  the  display  of  centerline  wanders.  The  centerline  was  shifted  from 
the  original  estimated  location  by  2.54  machine  centimeters  and  the 
total  integrated  power  recalculated  to  show  the  sensitivity  of  the 
centerline  location.  The  original  and  the  shifted  powers  differed  by 
+.01%.  The  largest  error  results  from  non-linear  and  oscillating  junc- 
tions between  digitized  data  points.  These  errors  are  estimated  to  be 
less  than  2.5%  in  incremental  area. 

Power  is  conserved  in  the  scattering  process.  If  the  digitiz- 
ing process  is  carried  out  sufficiently  far  for  both  the  turbulent  case 
as  well  as  the  non-turbulent  case,  the  total  power  is  the  same.  The 
runs  shown  in  Figures  42,  43,  and  44  have  power  discrepancies  of  -2.23% 
and  +10.14%,  and  -2.86%  in  total  power,  respectively,  for  Figures  42, 

43,  and  44  relative  to  the  power  shown  in  Figure  41.  Since  over  470 
traverses  were  reduced,  it  is  impractical  to  show  all  curves.  Figures 
41,  42,  43,  and  44  are  representative  of  the  least  well  behaved  traverse 
measurements.  In  general,  the  profile  measurements  become  better  be- 
haved as  L/Dp  increases  and  smaller  diameter  beams  are  employed,  as  may 
be  observed  by  comparing  Figures  39  and  40  with  Figures  41  and  43.  The 
power  discrepancies  would  also  decrease  with  beam  diameter  and  increase 
in  L/Dp. 

The  total  laser  power  was  normally  measured  without  the  aper- 
ture  in  front  of  the  power  meter  prior  to  collecting  a set  of  data.  The 
aperture  measured  and  calculated  powers  agree  with  the  total  laser  power 
measured  without  the  aperture  when  the  traverses  are  sufficiently  large. 
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In  digitizing  smooth  curves,  the  assumed  linear  fit  between  successive 
intensity  points  is  a very  good  approximation.  Some  turbulent  cases, 
with  large  values  of  scintillation,  resulted  in  oscillating  functions. 
Points  were  chosen  for  these  curves  such  that  differences  and  excesses 
were  approximately  equal  in  area. 
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CHAPTER  V. 


EXPERIMENTAL  RESULTS 

Summary 

In  this  chapter  the  mean  and  fluctuating  tunnel  profiles  will  be 
presented  for  velocity,  temperature,  index  of  refraction,  and  scales  of 
turbulence.  Next,  the  experimental  effects  on  a laser  beam  propagating 
through  this  defined  turbulent  medium  will  be  examined  by  observing: 

a.  The  distribution  of  beam  spread 

b.  The  increase  of  R/R  at  the  half  power  point  of  the 

0 

turbulent  laser  profile 

c.  The  effect  of  <An>2  on  laser  beam  spread 

d.  The  effect  of  X on  laser  beam  spread 

e.  The  effect  of  laser  beam  diameter  on  laser  beam  spread 

f.  The  effect  of  A on  laser  beam  spread 

g.  The  effect  of  diameter  of  the  laser  beam/A  ratio  on 

laser  beam  spread 

h.  The  reduction  of  centerline  intensity  of  the  laser  beam 
Table  2 illustrates  the  values  of  the  independently  variable  para 

meters  for  which  laser  beam  spread  measurements  were  obtained.  For 
example,  the  50  millimeter  diameter  helium  cadium  beam  spreao  was  in- 
ured at  L/0p  stations  of  28,  34.5  and  45  for  at  least  three  values  o,r 
heat  addition,  at  each  L/Dp station , for  each  tunnel  velocity  press"*-? 
at  which  significant  beam  spread  existed. 
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Table  2.  Range  of  Independent  Experimental  Parameters 
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Mean  and  Turbulent  Tunnel  Profiles 


Due  to  the  decreases,  primarily  in  <An>2,  which  results  in  reduced 
beam  spread,  the  number  of  spread  measurements  of  interest  decreases 
with  the  L/Dp  ratio.  The  combinations  of  interest  results  in  150  tunnel 
profiles  for  turbulence  generator  B which  were  reduced  in  accordance 
with  the  procedures  of  Chapter  IV.  A similar  set  of  spread  measurements 
were  performed  for  turbulence  generator  A at  an  L/ Dp  station  of  14. 

Sample  profiles  for  the  mean  and  fluctuating  components  of  U and  T, 
obtained  with  velocity  turbulence  generator  B,  are  shown  in  Figures  45, 
46,  47  and  48  for  a nominal  (tunnel  velocity)  pressure  of  .025  centi- 
meter of  water  and  heat  additions  of  zero,  516,  814  and  1196  watts, 
respecti vely.  The  rms  temperature  scales  of  Figures  46,  47  and  48  are 
10  times,  30  times  and  30  times,  respectively,  the  zero  heat  temperature 
scale  of  Figure  45.  The  mean  temperature  and  velocity  are  well  behaved 
with  little  structuring  of  the  mean  flow.  The  rms  profiles,  on  the 
other  hand,  show  rather  wide  variations  of  unpredictable  structuring 

i 

which  indicates  the  desirability  of  calculation  of  a spatial  mean  value 
for  the  rms  quantities.  Vertical  profiles  were  also  evaluated  for  turbu- 
lence generator  B at  three  L/Dp  stations,  three  heats  and  a tunnel 
velocity  pressure  of  .076  centimeter  of  water  to  ensure  that  the  free 
convective  effects  were  negligible.  Figure  49  shows  a vertical  profile 
obtained  for  a nominal  high  heat  (1166  watts)  at  a velocity  pressure  of 
.076  centimeter  of  water,  at  the  L/Dp  station  of  34.5,  indicating  negli- 
gible free  convective  effects  and  reasonably  uniform  rms  profiles  over 
the  central  portion  of  the  tunnel  through  which  the  beam  traverses. 
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raverse,  L/Dp  = 28,  High  Heat 


Figure  49.  Vertical  Tunnel  Traverse,  L /Dp  = 34.5,  High  Heat 


Variations  of  Fluid  Properties  Along  the  Tunnel  Axis 


Reduction  of  the  profiles  for  low  and  medium  heat  additions  at  the 
same  tunnel  velocity  pressure  of  .025  centimeter  of  water  results  in  a 
well  behaved  variation  with  L/ Dp ; in  temperature  fluctuation,  tempera- 
ture intensity,  and  index  of  refraction  turbulence  intensity  as  shown 
respectively  in  Figures  50,  51,  and  52.  The  upper  curve  in  Figure  50 
represents  the  nominal  high  heat  addition;  intermediate,  medium  heat; 
lower,  low  heat.  The  points  located  just  above  the  abscissa  for  these 
curves  represent  the  zero  heat  addition  case.  Since  n is  very  nearly  1, 
Figure  52  is  also  representative  of  <An>.  Figure  52  indicates  a range 
of  4 to  1 in  An  is  available  at  the  three  tunnel  stations  between  low 
and  high  heat  additions.  As  tunnel  velocity  is  increased,  AT  is  reduced 
allowing  the  range  in  <An>  between  zero  heat  addition  and  the  curves  of 
Figure  52  to  be  investigated,  giving  a total  available  range  of  1.4  x 
10'^  in  <An>2. 

The  velocity  turbulence  intensity  and  the  scales  of  turbulence  were 
found  to  be  less  well  behaved.  Temperature  scales  of  turbulence  will 
next  be  briefly  reviewed  and  the  experimentally  obtained  curves  of  ve- 
locity turbulence  intensity  and  the  temperature  integral  scale  of  turbu- 
lence given. 

Equation  (42)  gives  the  semi-empirical  relation  between  integral 
scale  of  turbulence  and  the  microscale  as: 


Batchelor  (Batchelor,  1953,  p.  136)  gives  the  decay  curve  for  t,Q2  for  the 
initial  period  as: 

1 0 v Dd 

V ■ < — D— E>  <86> 

Squaring  equation  (4<!)  and  substituting  the  iQ2  for  the  initial  period 
of  decay  gives  the  integral  scales  decay  as: 

10  Dp  x 

» = -f-Z  (o2)  (87) 

"g  P 

The  fluctuating  velocity  component  parallel  to  the  mean  velocity  is 
observed  to  decay  according  to  (Townsend,  1956,  p.  56): 


(u 1 )2  “ (x  " x0)n  (^8) 

where  n takes  on  the  values  from  1.2  to  1.0  for  non-isotropic  to  iso- 
tropic flows,  respectively,  or  for  isotropic  turbulences  in  terms  of  Tu 


1 U2  _ b ,x  ~ V ,aQ> 

TT^-^-  cD  < Dp  > m 

where  b is  a constant  dependent  upon  the  turbulence  generator  geometry 
and  Cp  is  the  drag  coefficient  of  the  generator.  b/Cp  may,  therefore, 
be  expected  to  be  weakly  Reynolds  number  dependent. 

The  prediction  of  the  integral  scale  and  the  turbulence  intensity's 
behavior  is  not  at  all  straightforward  as  is  illustrated  by  Batchelor's 
Figure  6.1  in  which  he  presents  the  graph  of  the  relation  (Batchelor, 
1953,  p.  106,  Fig.  6.1) 

dt  A 

or  du2  A , |qn 

dt-  Cu7  " (91) 
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This  relation  is  postulated  from  the  experimental  observation  that  the 
kinetic  energy  of  turbulence  decreases  as  u3/A  , with  C being  a constant 
of  order  unity.  Batchelor's  plot  of  this  relation  against  experimental 
observations,  as  a function  of  L/Dp  with  Reynolds  number  and  Dp  as  para- 
meters, shows  a spread  from  .8  to  1.45  in  du2A/(dt  Cu3)  at  an  L/Dp  of  20. 
Likewise,  there  is  not  a consistent  trend  with  Reynolds  number  change, 
nor  in  the  amount  of  change  in  du2A/(dt  Cu^with  Reynolds  number. 

The  behavior  of  the  scales  of  turbulence  may  be  expected  to  be 
extremely  difficult  to  predict  accurately  and  they  must  be  measured  for 
each  configuration  and  location  if  a value  is  to  be  obtained.  The  wide 
range  of  variation  of  the  integral  scale's  decay  is  additional  reason 
for  examination  of  the  microscale's  behavior  and  decay. 

The  zero  heat  addition  velocity  turbulence  intensity  is  shown  in 
Figure  53.  A spread  of  3 to  10%  in  the  values  of  Tu  with  Reynolds  num- 
ber is  observed  from  Figure  50.  The  .025  and  .076  centimeters  of  water 

velocity  pressure  Tu  points  are  well  behaved.  The  curves  through  these 

have  been  used  to  construct  the  mean  curve  for  the  decay  of  Tu  shown  in 

Figure  54  for  the  zero  heat  addition  case. 

The  velocity  Tu  generator  was  originally  envisioned  as  an  independ- 
ent control  on  the  integral  scale  of  the  temperature  fluctuations.  Com- 
parison of  the  temperature  integral  scales  with  the  velocity  integral 
scales  indicates  they  are  approximately  equal.  The  range  of  value  for 
the  velocity  scale  at  an  L/Dp  of  34.5  was  from  .953  to  1.499  centimeters. 
The  temperature  scales  at  the  34.5  station  ranged  from  .89  to  1.60  centi- 
meters. Comparison  of  the  velocity  and  temperature  integral  scales  at 
the  L/Dp  stations  of  28  and  45  were  also  similar  in  their  variations. 
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The  temperature  and  velocity  integral  scales  may  be  considered  equiva- 
lent in  magnitude.  The  temperature  integral  scale  is  shown  in  Figure 
54.  Examination  of  Figure  54  indicates  that  heat  addition  alters  the 
apparent  integral  scale  at  a given  L/Dp  location  in  an  irregular  fash- 
ion. The  trends  with  Reynolds  number  do,  however,  show  some  approximate 
ordering,  as  has  been  indicated  by  the  appropriate  curves,  for  velocity 
pressures  of  .025,  .051  and  .127  centimeters  of  water. 

Variations  of  Fluid  Properties  Along  the  Tunnel  Axis  with  the  Non- 
Dimensional  Heat  Addition  Parameter 

A method  for  normalizing  the  fluctuations  in  index  of  refraction  in 
the  tunnel  to  account  for  changes  in  tunnel  velocity  and  amounts  of  heat 
addition  which  varied  from  run  to  run  and  day  to  day  is  desirable,  par- 
ticularly at  low  tunnel  velocities  which  were  difficult  to  reproducibly 
set.  Equation  (27)  gave: 

^rms  , 

<An>  * -J2 — (92) 

while  T2  is  of  course  <*Q2,  T is  in  absolute  units,  and  for  this  experi- 
ment does  not  change  by  more  than  7%.  The  variation  in  <An>  will  be 
nearly  equivalent  to  the  change  in  dTrm$.  dTrms  is  bounded  by,  but  not 
equal  to  AT,  the  change  in  temperature  of  the  tunnel  flow.  The  assump- 
tion of  <An>  proportional  to  Q would  be  a reasonable  choice.  If  Q is 
non-dimensional ized,  the  parameter  becomes: 

t-9 (93) 

M Cp  T 

or  since  p remains  nearly  equal  to  barometric  pressure  in  this  experi- 
ment, as  was  pointed  out  in  Chapter  IV,  the  non-dimensional  parameter 


ill 


becomes: 


PB  Cp  A U " U 
Rn 


The  measured  values  of  ATrms»  T,  and  Pg  were  used  to  calculate  <An>  for 
the  integrated  tunnel  profiles  and  also  the  tunnel  centerline  values  in 
accordance  with  the  procedures  of  Chapter  IV.  The  mean  tunnel  traverse 
variation  of  <An>  with  the  non-dimensional  heat  addition  parameter  K Q/U 
is  shown  in  Figures  55,  56  and  57  for  turbulence  generator  B at  the  L/ Dp 
stations  of  28.  34.5  and  45  respectively.  A least  squares  linear  fit  of 
the  data  shows  the  non-dimensional  heat  addition  parameter  to  be  a 
reasonable  choice  for  comparing  and  predicting  the  <An>  variation.  At 
an  L/Dp  station  of  45,  Figure  58  shows  an  excellent  fit  with  K Q/U  for 
low,  medium,  and  high  values  of  heat  addition  at  approximately  the  same 
value  of  tunnel  velocity. 

Figure  59  shows  the  <An>  variation  with  K Q/U  for  turbulence  gener- 
ator A,  the  1.27  centimeter  perforation  diameter  turbulence  generator, 
at  an  L/Dp  station  of  14.  The  slope  of  the  <An>  curve  for  the  turbulence 
generator  A is  nearly  equal  to  turbulence  generator  B at  the  L/Dp  station 
of  28.  The  major  difference  between  generator  A and  B is  doubling  of 
the  turbulence  scale. 

The  velocity  turbulence  intensity  is  shown  in  Figures  60,  61  and 
62  against  the  parameter  K Q/U  at  the  L/Dp  stations  of  28,  34.5  and  45 
for  turbulence  generator  B.  The  velocity  turbulence  intensity  has  been 
calculated  to  include  the  temperature  correction  terms  as  given  by 
equation  (38)  of  Chapter  IV.  A linear  relation  is  again  shown  with  the 
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Figure  55,  </ir)>  versus  KQ/U  at  L/Dd  = 28,  Turbulence  Generator 


Figure  55.  <An>  versus  KQ/U  at  L/D 


.. 


Figure  57.  <An>  versus  KQ/U  at  L/Dp  = 45,  Turbulence  Generator 


Figure  58.  <An>  at  L/Dp  = 45  for  Low,  Medium  and  High  Heat 
Addition,  U = .99m/sec  to  1.16  m/sec 


<An>  versus  ICQ/U  at  L/'j0  = 14,  Turb 


Figure  60.  Centerline  Variation  of  Tu  at  L/D( 
Turbulence  Generator  B 


Centerline  Variation 
Turbulence  Generator 


non-dimensional  heat  addition  parameter  K Q/U.  The  temperature  correc- 
tion terms  are  very  important  in  determining  the  value  of  Tu.  The  Tu 


curve  for  the  turbulence  generator  A is  shown  in  Figure  63  at  an  L/Dp 
station  of  14.  At  the  L/ Dp  station  of  14,  considerable  structuring  of 
the  flow  is  observable  and  as  a result  comparably  more  data  spread  is 
present  in  Figure  63. 

The  scales  of  turbulence  were  previously  indicated  to  be  extremely 
difficult  to  predict.  The  measured  values  for  the  integral  scale  of 
temperature  are  shown  in  Figures  64,  65  and  66  against  K Q/U.  A general 
trend  may  be  seen  in  Figures  64,  65  and  66,  although  an  accurate  predic- 
tion is  not  possible.  Collapse  of  the  data  with  any  other  parameter 
would  not  be  expected  from  the  observations  of  this  data  and  also  from 
the  previously  indicated  variations  of  the  scales  (Batchelor,  1953, 
p.  106). 

The  rather  large  and  unpredictable  variation  of  the  integral  scale 
prompted  consideration  of  the  microscales'  behavior,  since  the  self- 
preserving character  of  the  correlation  function  should  allow  the  near- 
zero  time  behavior  to  be  constant  (Townsend,  1956,  p.  41).  The  micro- 
scales were  obtained  by  the  three  methods  indicated  in  Chapter  IV.  All 
of  the  methods  discussed  are  very  dependent  upon  the  specific  points 
chosen  and  the  character  of  the  correlation  function.  In  general,  the 

average  parabola  fit  gives  the  smallest  value  of  microscale,  the  fourth 

2 

order  polynominal  fit,  the  next  smallest  value,  and  the  C^  method,  the 
largest  value  of  microscale.  The  average  of  several  parabola  fits  was 
chosen  to  show  the  measured  microscales'  behavior  with  the  non-dimensional 
heat  addition  parameter,  since  in  absolute  value  it  shows  values  nearer 
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Figure  64.  Centerline  Value  of  the  Integral  Temperature  Seal* 
Turbulence,  L/C-,  = 28,  Turbulence  Generator  B 
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Lcnterhne  Value  ot  the  Integra!  Temperature 
Turbulence,  L/Dp  = 34.5,  Turbulence  Generator 


Figure  £6.  Centerline  Valup  of  the  Integral  Temperature  Scale  of 
Turbulence,  L/0n  = 45,  Turbulence  Generator  B 


the  empirical  predictions  of  equation  (85).  Figures  67,  68,  and  69 
show  the  behavior  of  the  microscale  with  K Q/U  at  the  L/Dp  stations  of 
28,  34.5,  and  45  respectively  for  turbulence  generator  B.  The  micro- 
scale Joes  appear  to  exhibit  more  regular  benav>or  at  tne  L/Dp  station 
of  28;  however,  tne  34.5  and  45  stafions  show  similar  ranges  and 
soreads  to  those  exhibited  by  the  integral  scale. 


Laser  beam  Spread  Measurements 

From  tne  axial  tunnel  profiles  just  presented,  the  prediction  of 
beam  spreac  with  <&r,>  may  readily  be  accomplished  for  the  general  range 
of  experimental  parameters,  while  the  deduction  of  the  effects  of  the 
scales  cf  turbulence  will  be  very  difficult  to  obtain  in  general  terms. 

Chernov's  calculation  of  the  focal  spot  size  will  be  used  to  assess 
the  predicted  behavior  at  the  focal  spot  in  the  presence  of  turbulence. 
Chernov  calculates  the  effective  spot  size  in  the  presence  of  turbulence 
to  the  non-turbulent  case  as  (Chernov,  1960,  p.  142): 


eff 


/a  h 


(95) 


Ay  7i  a 

where  u is  the  scattering  coefficient  given  by  (Sutton,  1969,  p.  1740): 

a = amplitude7  + phase7  (96) 


a = A <An>2  k2  a L0  (1  + j^f-~)'5/6 


(97) 


for  a gaussian  distribution  where  k«  is  <<1 , which  is  generally  the  case 
for  this  experiment,  a is  the  correlation  distance  or  scale  associated 
with  the  disturbance  - approximately  the  integral  scale.  The  value  of 
a for  a three-dimensional  turbulence  spectrum  is  given  by  (Sutton,  1969, 


L2b 


Turbulence,  L/Dp  = 34.5,  Turbulence  Generator 


69.  Centerline  Value  of  the  Micro  Temperature  Scale  of 
Turbulence,  L/Dp  = 45,  Turbulence  Generator  B 


p.  1740): 
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Substitution  of  « ir.-o  aqu:*v  i r '■ 


r - 2 * ■ n>  ! 

ert 
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task  trie  beam  1 


let*  r ■ . v>  f 


the  leaving  field  is  known.  . nov  speru  ,at--*s  i-; 


rea-' enable  tc 


assume  that  "the  field  near  the  ~er.  s o"  a ens  ii  ; ‘ in  its  icir 

features  with  the  field  of  a lens  1 :roi ' d b;  : - rrn  ar  diaper- :m 
Chernov's  square  aiaphrair.  was  required  to  c '.c  Vt  - • equired  ;”t2- 

arals  over  the  aperture.  In  this  worx  a rip  -i’  ir  aperture  wili  be  used 
and  the  mask  characteristic  area  hR?  will  be  replaced  b>  -u/4,  toe  mas< 
area,  and  the  area  of  the  beam  entering  the  turbulent  medium.  Chernov's 
expression  for  this  experiment  men  bacon's; 


3.54  <An>  D /L 


2 {(l+(f§£V  )'5/6}1/: 


(100) 


Sutton's  expression  for  tne  beam  spreao  at  the  half  power  point,  equa- 
tion (6),  is  the  same  as  equation  (100)  Since  nearly  !:  ::.i -t  ts  used 
in  laser  cavities  are  circular,  and  all  correspond ; -ig  a;  e -teres  are  cir 


cular,  the  evaluation  of  principal  r / ■„  dependencies  oi  a ‘ cam  r. p r 
tured  with  a circular  stop,  at  th c-  focus,  is  a’>so  very  important. 


Distribution  of  Beam  Spread 

The  emergence  of  a gaussian  beam,  as  a qaussian  beam  after  propa- 
gating through  a turbulent  medium,  is  not  a well  documented  assumption 


in  the  literature.  Therefore,  the  distribution  of  beam  spread  will  be 
examined  experimentally  to  determine  if  this  assumption  is  valid  and  at 
what  point  of  the  distribution  the  r^^/Ay  may  be  validly  applied,  if 
any.  Associated  with  the  distribution  of  intensity,  there  are  two 
problems  which  need  tc  be  addressed. 

The  first  arises  simply  in  comparing  runs  from  day  to  day,  for  geo- 
metric changes  in  beam  size,  ana  for  comparison  of  only  a portion  ot  tne 
total  distribution  to  a total  non-turbulent  distribution.  To  examine 
and  display  the  errors  which  are  introduced  by  truncation  of  the  distri- 
bution, where  the  traverse  is  not  carried  out  to  the  same  intensity,  a 
gaussian  distribution  was  digitized  in  accordance  with  the  procedures  of 

Chapter  IV.  The  resulting  curves  are  shown  in  Figure  70.  The  upper 

-x2 

curve  is  referenced  to  r (distance  from  Io)  at  values  of  e of  3.94  x 

10'4,  the  second  to  3 x 10  the  third  to  1.83  x 10  the  fourth  to 

7.7  x 10  , and  the  fifth  to  .61.  The  experimental  range  of  turbulent 

-x2 

to  non-turbulent  profiles  typically  ranged  from  .0183  to  .003  in  e 
The  error  introduced  by  not  obtaining  the  total  turbulent  distribution 
should  not  affect  the  runs  in  one  set  of  data  since  the  reference  beam 
remains  the  same.  However,  in  comparing  sets  of  data,  an  error  of  up  to 
approximately  10%  in  r at  the  half  power  point  may  be  introduced  when 

“X 

e at  the  tail  of  the  distribution  is  .003  instead  of  .000394. 

Because  each  set  of  runs  is  normalized  to  the  non-turbulent  case, 
valid  comparisons  are  possible  between  runs  and  any  errors  of  this  kind 
should  only  be  significant  if  there  has  been  a zero  shift  during  a run. 
Normally,  the  data  was  discarded  if  a significant  shift  in  zero  or  peak 
power  occurred. 


131 


r 


The  second  question  of  the  beam  intensity  distribution  is  whether 
the  beam  maintains  a gaussian  profile  after  traversing  the  turbulent 
medium,  as  is  assumed  in  many  analyses.  A limited  comparison  of  the 

I 

turbulent  profiles  to  a gaussian  profile  shows  differences  on  the  order 
of  10%,  as  shown  in  Figure  71,  in  I/Io.  The  upper  curve  in  Figure  71 
represents  the  turbulent  profile,  the  lower  the  gaussian.  Both  profiles 

I 

I 

were  integrated  to  the  same  I/Io  value  to  eliminate  the  effect  shown  in 
Figure  70.  Several  profiles  show  the  half  I/Io  points  and  1/n  I/Io 
points  remain  very  close  to  gaussian  values.  The  half  I/Io  points  were 
chosen  for  beam  spread  measurements.  The  measured  profile  was  shown  in 
Chapter  IV  to  be  dependent  on  the  aperture  which  is  selected  for  the 
traverse  of  the  beam  profile.  Any  conclusion  as  to  the  precise  varia- 
tion of  the  turbulent  profile  would  require  an  additional  investigation. 

Several  other  formulations  of  the  spot  size  have  been  published  in 
the  literature.  Sutton's  (Sutton,  1969,  p.  1741)  expression  for  the 
half  power  point  was  used  in  equation  (6),  to  illustrate  the  quantities 
required  in  evaluation  of  beam  spread.  The  functional  dependence  of 
Sutton's  equation  is  identical  to  Chernov's  for  the  case  when  D/A  >>1 . 

If  Sutton's  equation  is  parametrically  plotted,  a helium  cadmium  beam 
with  D/A  ranging  from  5.5  to  .5  in  increments  of  .05  in  D/A,  Figure  72 
is  obtained.  Figure  72  shows  that  at  the  half  power  point  of  the  dis- 
tribution there  is  no  discernible  curvature  introduced  by  the  -5/6  power 
term.  All  subsequent  plots  at  the  half  power  point  will,  therefore,  be 
shown  against  <An>  or  K Q/U. 

The  experimental  measurements  indicated  the  microscale  and  the  inte- 
gral scales  varied  with  heat  addition;  however,  there  were  cases  such  as 
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.076  centimeter  tunnel  velocity  pressure  under  which  the  integral  scale 
of  temperature  stayed  approximately  constant.  The  half  power  point 
variation  is  shown  in  Figures  73  and  74  for  two  tunnel  velocity  pres- 
sures. Figure  73  is  a least  squares  fit  through  the  zero,  low,  medium, 
and  high  heat  conditions  at  one  tunnel  velocity  pressure.  Figure  73  is 
for  a 50  millimeter  diameter  He  Cd  beam  at  the  L/Dp  station  of  28.  Fig- 
ure 74  shows  a similar  case  obtained  for  an  He  Cd  beam  of  30.58  milli- 
meters diameter  at  the  L/ Dp  station  of  45.  In  Figure  74  the  zero  heat 
point  at  .52  R/RQ  has  not  been  used  in  the  least  squares  fit  of  the  data 
shown  in  plotting  Figure  74.  The  zero  heat  point  could  normally  be  ex- 
pected to  remain  constant  within  +1%;  however,  this  variation  occasion- 
ally reached  +5%  due  to  thermal  storage  in  the  tunnel  walls  and  changes 
in  the  room  temperature.  Changes  in  room  temperature  can  alter  the 
laser  beam  power  entering  fluctuating  medium,  while  thermal  storage  in 
the  tunnel  walls  can  add  temperature  fluctuation  to  the  flow  just  as  the 
temperature  turbulence  generators  do.  The  effect  from  the  amount  of 
tunnel  heat  addition  was  normally  large  compared  to  the  effect  of  tunnel 
storage  and  room  temperature  changes.  In  general,  a data  set  for  a 
given  beam  diameter  and  L/Dp  station  would  yield  one  to  three  plots  such 
as  Figures  73  and  74  if  separately  plotted  for  each  tunnel  velocity 
pressure  as  heat  addition  is  varied. 

The  integral  scale  of  temperature  was  previously  indicated  to  be  a 
variable  parameter  when  heat  is  added  and  as  L/Dp  increases.  Equation 
(100)  shows  the  functional  dependence  of  beam  spread  is  a (/a)  depend- 
ence. The  plots  of  all  tunnel  velocity  pressure  operating  points  may  be 
expected  to  exhibit  some  data  spread  due  to  changes  in  A. 
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Figure  74.  Half  Power  Point  Variation  of  Focal  Spot  Size 
He  Cd,  L/DP  = 45,  D = 30.58  mm 


The  least  squares  fitted  curves  obtained  for  beam  spread  or  the  R/Rc 
ratio  at  the  .5  power  points  of  the  turbulent  and  non-turbulent  laser 
profiles  obtained  from  the  calculation  of  Chapter  IV  are  presented  in 
Figures  75  through  89.  At  the  .5  PQ  point  from  equation  (6) 


d (4X>  _ 3 MD  76,1/2 


' ay...:.  = 

d <Ao> 


° 

>.  /a  -75D 


Differentiation  of  the  two  linear  relations 


(101) 


=1 


d(jf)  = d (H  + J R/RQ) 


(102) 


d(Ao  x 10~7)  = d (F  + G {^) 

gives  upon  substitution  of  equation  (102)  into  equation  (103) 

d(R/R0)  1 

d <An>  GJ 

which  also  may  be  equated  to  equation  (101)  giving: 

r.  1./2ttA»2,-5/6-i1/2  D Ap  _ _J 

L ;75D  J , rr  3.54GJ 


(103) 


(104) 


(105) 


Equation  (105)  may  be  used  with  the  mean  profile  data  at  the  L/Dp  sta- 
tion of  28  for  the  He  Cd  and  He  Me  data  to  compare  the  wave  length  de- 
pendence of  equation  (105). 

Taking  the  ratio  of  equation  (105)  for  He  Cd  to  He  Ne: 


{i+(^)2V5/6T/z 

L'  .750  ’ J 

HHfcT5/6]1/s 


1/2 

\X  /a  He  Cd  _ (G  J)  He  Ne 
1/2  / D/Ln\  ^ He  Cd 


A /A  / He  Ne 


(106) 


or  since  A/D  is  approximately  the  same  for  the  He  Ne  and  He  Cd  beams: 

(107) 

The  values  experimentally  obtained  are  given  in  Table  3. 


(x)  He  Ne  _ (J)  He  Ne 
(xTHe~Cd  (J)  He  Cd 


Table  3.  x Dependence  on  Beam  Spread 


Diameter  of 
Beam-Laser 

Experimental 
Value  of 
A He  Ne/xHe  Cd 

Calculated 
Value  of 
xHe  Ne/xHe  Cd 

50mm-He  Ne 
50mm-He  Cd 

1 .375 

1.433 

26mm-He  Ne 
30.58mm-He  Cd 

1 .660 

1 .685 

8mm-He  Ne 
14.03-He  Cd 

2.660 

2.513 

Next,  if  the  ratio  of  equation  (105)  is  taken  for  the  beam  diameter  at 
each  station,  the  apparent  a scale  and  x are  constant,  as  is  the  value 
of  the  slope  of  the  <An>  versus  K Q/U  curve,  G.  The  ratio  of  J's  or  the 
slopes  of  the  R/RQ  curves  at  .5  PQ  at  the  same  wave  length  then  gives 
the  beam  diameter  ratio  dependence.  Taking  the  ratio  of  equation  (104) 
to  the  same  equation  for  the  50  millimeter  diameter  beam  gives: 

in  * OV5/Y/2  „ 


: 75D; 

r 2tt  A 


50 


in  - (fit)  > ' I 


-5/6. 1/2  D 


50 


(108) 


50 


Neglect  of  a/D  terms  results  in  the  experimental  and  actual  values  which 
are  presented  in  Table  4.  The  *value  was  computed  for  a set  of  <An>  vs. 

K Q/U  data  which  had  a very  linear  experimental  range  and  the  correspond- 
ing R/Rq  data  also  had  a corresponding  linear  range  as  is  shown  in  Fig- 
ures 90  and  91 . 
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Figure  75.  Half  Power  Point  Variation  of  Focal  Spot  Size, 
He  Cd,  L/Dp  = 28,  D = 50  mm 


Figure  79.  Half  Power  Point  Variation  of  Focal  Spot  Size, 
He  Cd,  L/Dd  = 34.5,  D = 50  mm 


igure  80.  Half  Power  Point  Variation  of  Focal  Spo 
He  Cd,  L/Dd  = 34.5,  D = 30.58  mm 


Half  Power  Point  Variation  of  Focal  Spot  Size 
He  Cd,  L/Du  = 34.5,  D * 14.03  mm 


Figure  82.  Half  Power  Point  Variation  of  Focal  Spot  Size 
He  Cd,  L/Dd  = 45,  D = 50  mm 


Figure  83.  Half  Power  Point  Variation  of  Focal  Spot  Size, 
He  Cd,  L/Dd  = 45,  D = 30.53  mm 


Figure  87.  Half  Power  Point  Variation  of  Focal  Spot  Size, 
He  Ne,  L/D„  = 28,  D = 26  mm 


Table  4.  Diameter  Dependence  of  Beam  Spread 


X 

L/Dp 

Calculated 
Value  J50/J50 

Actual 

Value 

He  Cd 

28 

1.000 

1 .000 

441 .6mm 

.920 

.616 

.451 

.281 

34.5 

1.000 

1.00C 

.817 

.616 

.363 

.281 

45 

1.000 

1.000 

.697* 

.616 

.823 

.616 

.157 

.281 

He  Ne 

28 

1.000 

1 .000 

632.8mm 

.766 

.520 

.234 

.160 

*Value 


The  diameter  dependence  is  not  as  good  as  the  wave  length  depend- 
ence because  there  occur  ranges  of  experimental  data  for  which  D/A  is  of 
order  1 and  the  neglected  term  is  important.  Under  these  conditions, 
there  is  little  or  no  beam  spread  since  the  beam  sees  nearly  uniform 
conditions  as  the  distorted  medium  passes,  with  only  centroid  motion  and 
scintillation  resulting.  For  the  beam  14.03  millimeters  diameter,  the 
turbulence  scale  approaches  beam  size  and  a resulting  near  vertical  por- 
tion of  the  R/Rq  at  .5  PQ  versus  K Q/U  occurs,  distorting  the  slope  at  a 
given  value  of  An.  The  curvature  of  the  data  would  indicate  that  a sim- 
ple D variation  is  not  correct. 

From  equation  (95)  the  effective  scale  for  a gaussian  distribution 
may  be  calculated  which  would  produce  the  measured  beam  spread.  Corre- 
lation functions  were  obtained  in  this  experiment  which  were  typical  of 
a single  scale,  three  different  scales,  and  even  functions  typical  of  an 


infinite  range  of  scales.  A single  scale  is,  therefore,  not  a realistic 
representation  of  the  flow  field  actually  encountered  in  this  experiment. 
The  implied  effective  gaussian  scales  are  tabulated  in  Table  5 for  the 
mean  values  of  the  tunnel  profiles  previously  presented. 

Table  5.  Calculated  Gaussian  Scale  from  Mean  Tunnel  Profiles 


X 

L/Dp 

Beam  Diameter 

a 

He  Cd 

28 

50.00mm 

29.7mm 

441 .6mm 

30.58 

13.0 

Turbulence 

14.03 

11.6 

Generator  B 

1.50 

No  Spread 

34.5 

50.00 

3.0 

30.58 

2.5 

14.03 

1 .1 

45 

50.00 

2.8 

30.58 

1.9 

14.03 

2.3 

1.50 

No  Spread 

He  Ne 

28 

50.00 

3.1 

632.8mm 

26.00 

1.1 

Turbulence 

8.00 

.7 

Generator  B 

He  Ne 
632 . 8mm 
Turbulence 
Generator  A 

14 

50.00 

3.0 

The  single  gaussian  scale  is  observed  to  be  in  the  range  of  the 
measured  integral  and  microscale  shown  previously.  The  very  small  value 
of  (a)  indicates  very  little  beam  spread.  Integral  scales  from  6.4  to 
2.3  millimeters  were  measured  for  turbulence  generator  B. 

Turbulence  generator  A had  twice  the  perforation  diameter  of  turbu- 
lence generator  B which  resulted  in  the  integral  scale  being  comparable 
to  the  beam  diameter  for  the  50  millimeter  diameter  beam.  Very  little 
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spread  would  be  predicted  for  the  50  millimeter  beam.  The  R/RQ  spread 
at  .5  P0  point  for  the  He  Ne  beam  and  turbulence  generator  A is  shown  in 
Figure  92,  further  illustrating  the  difficulty  in  using  a single  effec- 
tive scale  when  the  integral  scale  approaches  the  beam  diameter. 

As  a result  of  the  inadequacies  shown  in  Tables  4 and  5,  which  occur 
with  neglect  of  the  A/D  terms  of  equation  (100),  the  experimental  laser 
beam  spread  at  the  half  power  point  has  been  shown  in  Figures  93  through 
105,  along  with  the  predicted  values  of  spread  for  equation  (100)  as  Ay 
is  varied.  The  solid  curves  in  Figures  94  through  104  are  from  equation 
(100),  with  the  integral  scale,  Ay,  varied  from  .254  cm  in  .254  cm  incre- 
ments. The  curve  nearest  the  ordinate  represents  equation  (100)  for  an 
integral  scale  of  .254  cm.  In  Figure  105,  equation  (100)  is  shown  for 
Ay  of  .127  and  .254  cm.  Quantitatively,  the  diameter  dependence  shown 
in  Figures  94  through  105  is  much  better  than  indicated  in  Table  4.  The 
integral  scale  dependence  is  also  shown  to  be  in  the  range  of  the  experi- 
mentally measured  values.  Equation  (100)  does  appear  to  overestimate 
the  integral  scale  and  not  be  capable  of  accurately  predicting  the  spread 
for  very  small  values  of  <An>-  The  larger  values  of  < An>  tend  to  have 
smaller  measured  values  of  the  integral  scale  of  turbulence.  Figures  94 
through  105  also  show  an  increase  in  slope  with  <An>  or  indicate  a reduc- 
tion in  the  integral  scale  with  increasing  <An>.  Figure  105  indicates  a 
turbulence  scale  of  ^.178  cm.  The  microscale  of  turbulence  was  indicated 
to  range  from  .12  to  .5  cm  in  Figures  67  to  69.  The  values  of  Ay,  which 
are  implied  when  the  experimental  data  for  beam  spread  at  the  half  power 
point  of  the  power  distribution  are  compared  against  equation  (100),  are 
summarized  in  Table  6. 
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Figure  92.  Half  Power  Point  Variation  of  Focal  Spot  Size, 
He  Ne,  L/’Dd  = 14,  D = 50  mm 


Figure  97.  Half  Power  Point  Variation  for  the  Chernov 
Calculation,  He  Cd,  L/Dp  = 34.5,  D = 30.58 
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Figure  ICO.  Half  Power  Point  Variation  for  the  Chernov-Sutton 
Calculation,  He  Cd,  L/Dp,  = 45,  D = 30.58  mm 
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Figure  101.  Half  Power  Point  Variation  for  the  Chernov-Suttor 
Calculation,  He  Cd,  L/Dp  = 45,  D = 14. C3  mm 


Figure  102.  Half  Power  Point  Variation  for  the  Chernov-Sutton 
Calculation,  He  Ne,  L/Dd  = 34.5,  D = 50  nn 


Figure  103.  Half  Power  Point  Variation  for  the  Cnernov-Sutton 
Calculation,  He  Ne,  L/Dd  = 34.5,  D = 26  mm 


Table  6.  Comparison  of  Ay  and  i0  Scales  from  Half  Power  Values 
of  R/Ro  Measurements  with  Their  Measured  Values 


The  Reduction  of  Centerline  Intensity 


The  turbulent  to  non-turbulent  beam  intensity  ratio  is  given  by 
(Chernov,  1967,  p.  142): 

= tt  a2/ah02  (109) 

or  for  this  experiment  using  Sutton's  equation  (6)  for  a and  replacing 


a2  and  hQ2  with  (A/.75)2,  the  value  for  the  three-dimensional  turbulence 
spectrum  (Sutton,  1969,  p.  1740),  and  irD2/4  respectively: 


I 

Io 


= 

2(  .75  DTT<An>)2  Lo  (1  +(^-)2)"5/6 


(110) 


This  relation  is  plotted  parametrically  against  A/D  in  Figure  105  for  a 
50  mm  diameter  He  Cd  beam.  A/D  takes  on  values  from  .025  to  .5,  with 
the  .025  curve  being  closest  to  the  abscissa.  Figure  105  and  equation 
(110)  show  that  the  expression  1/ Io  does  not  go  smoothly  to  1.0  as  <An> 
goes  to  zero. 

The  experimentally  measured  centerline  intensity  ratios  are  shown 
in  Figures  106  to  117.  The  solid  curves  are  the  predicted  values  from 
equation  (110),  for  integral  scales  of  .254  cm  in  .254  cm  increments  for 
Figures  106  to  116.  The  solid  curves  of  Figure  117  begin  with  i.Q  of  .08 
cm  with  increments  of  .0254  cm  between  successive  curves.  Comparison  of 
the  similar  diameter  beams  at  the  same  <An>  values,  as  L/Dp  is  increased, 
shows  the  increase  in  the  integral  scale  of  turbulence  which  was  shown 
in  Figure  54.  The  integral  scale  of  turbulence  which  is  implied  in  com- 
parison of  the  experimental  data  with  the  parametric  plot  of  equation 
(110)  shown  in  Figures  106  through  117  has  been  tabulated  in  Table  6 for 
comparison  against  the  measured  temperature  scales  of  turbulence  from 
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Value  of  Intensity  for  a He  Cd  Beam 


Figure  107.  Centerline  Value  of  Intensity  for  a He  Cd  Beam 


Figure  108.  Centerline  Value  of  Intensity  for  a He  Cd  Beam, 
• L/Dd  = 28,  D = 14.03  rm 


Figure  109.  Centerline  Value  of  Intensity  for  a He  Cd  Beam 
L/Dn  = 34.5,  D = 50  mm 


Figure  110.  Centerline  Value  of  Intensity  for  a He  Cd  Beam 
L/Dn  = 34.5,  D = 30.53  mm 


Centerline  Value  of 
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Figure  115.  Centerline  Value  of  Intensity  for  a He  Ne  Beam 
L/Dd  = 34.5,  D = 20  mm 


Figure  117.  Centerline  Value  of  Intensity  for  a He  Ne  Beam 


Figures  64  through  69.  Figures  64  to  66  show  that  the  smaller  value  of 
Aj  given  in  the  table  occurs  for  the  larger  values  of  K Q/U  or  <An>. 

The  50mm  diameter  beams  show  very  good  agreement  with  the  predicted 
integral  scales  of  turbulence.  Comparison  of  the  50mm  He  Ne  data  with 
the  50mm  He  Cd  indicates  excellent  agreement  with  the  wave  length  de- 
pendence of  equation  (110).  As  the  ratio  of  beam  diameter  to  turbulence 
scale  approaches  1,  the  effective  turbulence  scale  approaches  the  micro- 
scale of  turbulence.  The  range  of  measured  values  of  the  temperature 
microscale  are  given  in  the  last  column  of  Table  7.  The  He  Cd  14.03mm 
diameter  beam  at  L/Dp  - 28  and  the  8mm  diameter  He  Ne  beams  show  values 
of  turbulence  scales  very  close  to  the  minimum  measured  microscale  and 
the  L/Dp  stations  of  34.5  and  45  for  the  14.03mm  diameter  He  Cd  beams 
indicate  microscales  near  the  maximum  indicated  values  measured  values. 
Replacement  of  the  integral  scale  with  the  microscale,  as  D/A  approaches 
1,  appears  approximately  correct.  Sutton's  expression  for  a is  shown  to 
give  a good  functional  description  of  I/Io  after  I/Io  is  below  75  to 
80%.  Quantitatively,  the  predicted  values  are  quite  good  below  the 
I/Io ' s of  .75  to  .80%. 
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CHAPTER  VI 


CONCLUSIONS 

The  objective  of  this  dissertation  was  to  determine  whether  the 
discrepancies  between  analytic  and  experimental  results  were  caused  by 
inaccuracies  of  the  analytic  models  or  by  difficulties  encountered  in 
making  precise  experimental  measurements  for  ranges  of  geometric  condi- 
tions typical  of  laser  cavities.  The  objective  required  the  creation 
of  controlled  turbulence  conditions  which  would  produce  beam  spreads  on 
the  order  of  1.5  to  2 times  that  of  a diffraction  limited  beam  diameter. 
In  this  dissertation,  experimental  beam  spread  measurements  have  been 
reported  for  the  geometric  range  of  conditions  encountered  in  and  around 
laser  cavities  with  the  ratio  of  the  integral  temperature  scale  of 
turbulence  to  the  laser  beam  diameter,  Aj/D,  varied  from  <1  through  >1. 
The  values  of  beam  spread  have  been  presented  for  the  range  of  turbu- 
lence Reynolds  numbers  which  are  normally  encountered  in  laser  cavities. 
The  attenuation  or  extinction  coefficient,  cxL0,  was  varied  from  <<1  to 
>1.  The  parameters  required  to  predict  laser  beam  spread  are  wave 
length,  laser  beam  diameter,  integral  turbulence  scale,  micro  turbulence 
scale,  and  the  fluctuation  of  the  index  of  refraction.  The  beam  spread 
for  variation  of  the  important  parameters  was  reported  at  the  half  power 
points  and  the  centerline  of  the  laser  intensity  distribution  for  wave 
lengths  of  441. 6nm  and  632. 8nm. 
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Crv  tov's  beam  spread  and  centerline  intensity  expressions  for  a 
rectangular  aperture  were  used  for  comparison  with  the  experimental  data. 
Chernov's  expression  was  found  to  apply  to  a circular  aperture  when  the 
beam  area,  irD2/4,  replaces  the  rectangular  aperture  area,  hQ2,  and  the 
three-dimensional  integral  scale,  a^./.75,  is  used  instead  of  the  energy 
scale,  a. 

The  laser  intensity  distribution  in  a turbulent  flow  field  was 
shown  to  be  approximately  gaussian.  The  normalized  value  of  the  focal 
spot  radius,  R/Ro,  at  the  half  power  and  1/e  points  of  the  power  distri- 
bution curve  is  normally  within  5 to  10%  of  the  gaussian  distribution 
when  all  distributions  are  normalized  to  the  same  normalized  intensity 
value,  I/Io. 

The  relationship  between  the  fluctuation  in  index  of  refraction  and 
the  value  of  R/Ro  at  the  half  power  point  with  the  rms  fluctuation  in 
index  of  refraction,  <An>,  was  found  to  be  linear,  within  the  measure- 
ment accuracy  of  individual  parameters,  and  significantly  better  than 
previously  reported  in  the  literature.  When  the  integral  scale  of  turbu- 
lence, Ay,  remained  constant,  deviations  from  linearity  were  within  +1/2%. 
Although  the  extinction  coefficient,  ciLq,  ranged  up  to  4 for  the  reported 
data,  the  expression  for  aLQ  <1  was  in  excellent  agreement  with  the  ex- 
perimental data  for  all  otLQ ' s with  values  of  I/Io  <.8. 

The  wave  length  dependence  of  beam  spread  was  shown  by  taking  the 
ratio  of  the  values  of  R/Ro  for  the  half  power  point  beams  of  helium  neon 
and  helium  cadmium.  The  experimental  beam  wave  length  ratio  of  helium 
neon  to  helium  cadmium.  AHe  Ne/AHe  Cd,  from  the  beam  spread  measurements 
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was  compared  in  Table  3 to  the  actual  wave  length  ratio.  Table  3 indi- 
cated agreement  of  the  experimental  and  actual  values  within  -4%,  -1%, 
+6%. 

The  beam  diameter  dependence  of  beam  spread  at  the  half  power  point 
obtained  for  Chernov's  equation  (95)  is  qualitatively  correct,  but  does 
not  predict  the  beam  diameter  effect  on  beam  spread.  Equation  (95), 
when  based  upon  beam  spread  measurements,  tends  to  overestimate  the  di- 
ameter dependence  by  a +50%  with  the  best  value  showing  a +13%  differ- 
ence when  a single  scale  of  turbulence  model  is  used. 

When  Chernov's  equation  (95)  is  modified  by  the  substitution  of  the 
laser  beam  area,  the  a,  and  Ay  values  for  a three-dimensional  spectrum 
of  turbulence,  a significant  improvement  in  the  quantitative  prediction 
of  laser  beam  spread  is  obtained.  The  best  previous  prediction  of  beam 
spread  compared  to  theory  has  been  indicated  to  only  be  within  a factor 
of  2 to  10.  Given  a value  of  <An>»  the  factor  of  2 uncertainty  in  the 
normalized  R/Ro  value  results  in  an  uncertainty  in  Ay  of  the  order  of 
+2 Ay , —1/2 Ay • The  largest  differences  from  the  predicted  R/Ro  values, 
which  were  observed  at  the  half  power  values  of  R/Ro,  were  -20%  to  -30% 
for  all  beam  diameters  except  the  8mm  helium  neon  beam  for  which  Ay/D 
ranges  from  1 to  2.  In  this  case,  the  value  of  the  turbulence  scale 
from  the  spread  measurements  is  the  same  as  the  measured  value  of  micro- 
scale, as  was  shown  in  Table  6.  Table  6 also  illustrates  the  excellent 
agreement  by  comparing  the  calculated  Ay  value  obtained  from  beam  spread 
measurements  with  the  measured  value  of  Ay. 

Examination  of  the  centerline  intensity  ratio,  by  means  of  Sutton's 
modified  a in  Chernov's  expression,  resulted  in  very  good  predictions  of 
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Ay/D  <.01  and  I/Io  <80%.  The  <Ar)>2,  Ay/D,  and  A trends  of  beam  spread 
were  very  good.  The  comparison  of  the  experimental  centerline  intensi- 
ties with  parametric  plots  of  equation  (110)  showed  excellent  agreement 
with  <An>2.  The  comparison  of  the  50mm  diameter  He  Ne  and  He  Cd  beam 
spreads  at  the  L/Dp  station  of  34.5  showed  excellent  agreement  for  the 
wave  length  dependence  of  equation  (110).  Sutton's  expression  for  the 
three-dimensional  spectrum  of  turbulence  gave  excellent  quantitative 
agreement  with  the  actual  measured  integral  scales  of  turbulence  for  the 
50mm  diameter  beams.  Integral  scales  of  turbulence  were  derived  from 
the  beam  spread  measurements  which  were  within  2 to  6%  of  the  actual 
values  at  all  stations  except  at  L/Dp  = 45.  Experimental  values  shown 
in  Figures  94  through  105  indicate  excellent  agreement  with  analytic 
values  for  these  conditions.  The  quantitative  agreement  for  Ay/D  <.01 
for  centerline  values  was  shown  to  be  very  good  and  better  than  any  data 
reported  previously. 

As  Ay/D  approache'  2,  it  appears  that  the  temperature  microscale 
of  turbulence  may  be  useo  to  replace  the  integral  scale.  The  value  of 
£Qy  obtained  from  beam  spread  measurements,  as  shown  in  Table  6,  is 
within  11%  of  the  measured  value  for  the  8mm  He  Ne  beam. 

For  intermediate  values  of  Ay/D  (approaching  .3),  equation  (110) 
is  found  to  be  inadequate  for  the  prediction  of  beam  spreads.  The  dis- 
crepancy in  temperature  scale  at  Ay/D  of  .2  is  20%,  and  34%  at  a Ay/D 
of  .3  for  the  He  Cd  data. 

Comparison  of  Tables  6 and  7,  which  summarizes  the  beam  spread  cal- 
culated and  measured  value  of  Ay,  indicates  that  measurements  of  laser 
beam  spread  at  the  half  power  point  «f  the  laser  power  distribution  are 
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more  accurate  in  describing  the  A^/D  variation  of  laser  beam  spread  by 
a ratio  of  nearly  2 to  1.  The  measurement  of  beam  spread,  for  all 
values  of  Ay/D  <1  at  the  half  power  point,  has  been  shown  to  agree  with 
the  predictions  of  equation  (100). 

The  results  of  this  dissertation  have  reduced  the  uncertainty  in 
predicting  laser  beam  spread  by  a factor  greater  than  2 for  a wide  range 
of  turbulence  conditions  which  are  encountered  in  laser  cavity  flows. 

As  a result  of  this  improvement  of  the  relationship  between  experiment 
and  theory,  the  uncertainty  in  predicting  losses  and  laser  beam  spread 
has  likewise  been  reduced  for  all  mediums  and  conditions  for  which  dif- 
fraction effects  may  be  neglected.  The  combination  of  Chernov's  and 
Sutton’s  relations  was  shown  to  accurately  predict  the  centerline  inten- 

J 

sity  variation  and  accurately  predict  the  half  power  point  variation  of 
beam  spread. 

The  primary  difficulty  in  predicting  beam  spread  caused  by  turbu- 
lence, in  the  range  of  parameters  investigated,  consists  in  accurately 
measuring  the  turbulence  parameters.  As  a result,  the  2 to  10  fold 
discrepancies,  which  appear  in  the  literature  in  predicting  the  beam 
focal  spot  size  resulting  from  turbulence,  may  be  attributed  primarily 
to  a lack  of  precise  measurements. 
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CHAPTER  VII 


RECOMMENDATIONS 

In  and  around  laser  cavities,  the  turbulent  environment  is  much  as 
was  simulated  in  this  experiment.  A/D  is  of  order  1 and  is  not  neg- 
ligible compared  to  A.  These  conditions  violate  several  of  the  analyti- 
cal solutions  available  in  the  literature.  It  is  recommended  that  the 
experimental  data  be  compared  quantitatively  against  available  numerical 
solutions  such  as  Lutomirski's  (Lutomirski , 1970)  or  Sutton's,  D/A  ^1, 
(Sutton,  1969).  A comparison  of  similar  general  conditions  indicates 
that  the  specific  agreement  may  be  quite  good. 

It  is  recommended  that  a similar  experiment  be  carried  out  for  a 
longer  wave  length,  such  as  CO^,  under  well  controlled  and  measured 
turbulence  conditions. 
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more  accurate  in  describing  the  A^/D  variation  of  laser  beam  spread  by 
a ratio  of  nearly  2 to  1 . The  measurement  of  beam  spread,  for  all 
values  of  A-J-/D  <1  at  the  half  power  point,  has  been  shown  to  agree  with 
the  predictions  of  equation  (100). 

The  results  of  this  dissertation  have  reduced  the  uncertainty  in 
predicting  laser  beam  spread  by  a factor  greater  than  2 for  a wide  range 
of  turbulence  conditions  which  are  encountered  in  laser  cavity  flows. 

As  a result  of  this  improvement  of  the  relationship  between  experiment 
and  theory,  the  uncertainty  in  predicting  losses  and  laser  beam  spread 
has  likewise  been  reduced  for  all  mediums  and  conditions  for  which  dif- 
fraction effects  may  be  neglected.  The  combination  of  Chernov's  and 
Sutton's  relations  was  shown  to  accurately  predict  the  centerline  inten- 
sity variation  and  accurately  predict  the  half  power  point  variation  of 
beam  spread. 

The  primary  difficulty  in  predicting  beam  spread  caused  by  turbu- 
lence, in  the  range  of  parameters  investigated,  consists  in  accurately 
measuring  the  turbulence  parameters.  As  a result,  the  2 to  10  fold 
discrepancies,  which  appear  in  the  literature  in  predicting  the  beam 
focal  spot  size  resulting  from  turbulence,  may  be  attributed  primarily 
to  a lack  of  precise  measurements. 
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